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SYNOPSIS 

This  report  presents  two  possible  methods  of  structural 
analysis  for  long  span  cable-stiffened  structures. 

The  first  method  is  to  analyze  the  cable  stiffened 
structure  by  the  principle  of  superposition,  in  which  the 
first  step  is  to  reduce  the  structure  to  be  a  statical  deter- 
minate or  statical  indeterminate  in  minimum  degrees  base 
structure  and  put  all  the  redundants  as  external  forces  to 
be  solved.   In  this  analysis,  redundants  are  the  cable  stresses. 
The  second  step  is  to  find  the  deflections  of  the  base  struc- 
ture at  each  cable  attachment  due  to  actual  loads  and  unit 
loads  by  the  method  of  moment  area  integration,  based  on 
Castigliano' s  theorem  that  the  deflection  at  the  cable  at- 
tachment is  the  summation  of  deflection  due  to  actual  loads 
and  the  unit  deflections  caused  by  unit  stress  on  cables  times 
the  cable  stresses.   Therefore,  a  number  of  simultaneous 
equations  can  be  found  by  the  equilibrium  of  cable  stress  and 
deflection  at  each  cable  attachment.   In  order  to  solve  those 
linear  equations,  the  author  has  used  matrix  method  analysis 
and  set  up  computer  programs  for  those  moment  area  integration, 
deflections  and  matrix  inversions  to  find  cable  stresses  in 
a  quick  way. 

The  second  method  of  analysis  is  to  assume  that  the  cables 
work  as  part  of  the  frame  reinforcement  and,  using  the  prin- 
ciple of  transformed  areas,  the  frame  can  be  interporated  as 
a  haunched  member  with  variable  moment  of  inertias.   That  is, 


the  section  modules  at  each  end  of  the  frame  is  increased 
due  to  cable  stiffening  without  increasing  concrete  quantities. 
Therefore,  the  cable  stiffened  structure  can  be  analyzed  by- 
using  general  haunched  member  theory.   The  author  used  a  com- 
puter and  solved  this  time  consuming  calculation  problem 
within  five  minutes. 

Numerical  examples  are  given  to  illustrate  those  methods 
and  several  computer  FORGO  programs  are  written  to  solve  those 
examples  with  high  accuracy.   Another  new  idea  developed  by 
the  author  is  using  prestress  for  a  further  reduction  of  end 
mement  of  the  frame.   That  is,  the  cable  can  be  pretensioned 
in  a  definite  amount  for  a  definite  amount  of  mement  to  be 
reduced  from  the  frame.   This  method  has  been  shown  in 
numerical  example  1,  which  is  helpful  in  actual  design. 


INTRODUCTION 

In  any  k.ind  of  structural  design,  there  are  two  essential 
conditions  for  economy  to  be  considered:   the  efficient  utili- 
zation of  material,  and  the  maximum  reduction  of  the  dead 
weight  of  the  structure. 

For  a  long  spanned  rigid  frame,  the  very  large  corner 
stresses  at  the  fixed-end  require  that  this  section  be  ex- 
tremely large.   This  means  not  only  that  the  concrete  and 
reinforcing  steel  shall  be  increased  but  also  the  dead  loads 
of  the  frame.   If  the  span  of  the  rigid  frame  is  very  long, 
it  is  almost  impossible  to  construct  without  providing  addi- 
tional support  between  the  span  to  reduce  the  fixed  end 
sectional  areas.   For  an  economical  design,  it  is  possible  to 
use  cable  stiffening,  because  a  cable-stiffened  structure 
not  only  can  make  a  very  large  space  without  intermediate 
support  but  also  can  reduce  the  dead  loads  as  well  as  amounts 
of  construction  materials. 

The  first  developed  cable  stiffened  structures  were  for 
the  girder  bridge  stiffening.   In  the  time  of  World  War  II, 
there  was  an  extreme  shortage  of  construction  materials  to 
rebuild  the  German  bridges  which  had  been  destroyed  in  the 
war.   So  German  engineers  designed  a  new  type  of  construction 
which  would  be  safe,  durable,  economical  in  cost  and  in  quan- 
tity of  construction  materials,  and  pleasing  architecturally 
all  at  the  same  time.   The  application  of  cable  stiffening 
to  architectural  structures  has  been  developed  only  in  the 


last  ten  or  twelve  years.   The  typical  example  of  cable 
stiffened  building  in  the  United  States  is  the  Pan  American 
Terminal  Building  at  Kennedy  International  Airport.   The 
purpose  underlying  the  design  is  to  permit  aircraft  to  come 
underneath,  for  maintenance  in  the  case  of  the  hangar,  and 
for  embarkation  of  passengers  in  the  case  of  this  terminal 
building.   In  both  schemes,  cables  perform  an  essential  role 
in  enabling  the  tremendous  cantilivers  to  be  built.   The 
silhouette  is  like  a  vast  umbrvella  covering  the  activities 
of  travel.   Many  other  buildings  have  been  constructed  in 
recent  years  in  cable  stiffening  forms  or  cable  suspended 
forms  to  make  the  structures  in  good  architectural  shape  and 
economical  in  material  utilization,  which  is  almost  impossible 
in  ordinary  construction  methods. 

In  this  report,  the  author  has  taken  a  typical  long  span 
rigid  frame  with  cable  stiffening  as  an  example  to  show  how 
to  analyze  the  cable  stiffened  structures.   Many  possible 
methods  of  analysis  can  be  employed  in  analyzing  this  type  of 
structure.   Here  the  author  introduces  two  kinds  of  methods 
of  structural  analysis  which  have  not  been  used  previously  in 
the  cable  stiffened  structural  analysis. 

The  first  one  is  to  solve  the  cable  stiffened  structure 
by  means  of  the  superposition  principle,  which  is  generally 
known  as  the  linear  structural  analysis  and  is  used  for  solving 
high  order  statically  indeterminate  structures.   In  this  study, 
the  cable  forces  are  taken  as  the  indeterminate  redundants  to 


be  solved  with  the  frame  treated  as  a  base  structure  without 
cable  stiffening.   The  moment  on  the  base  structure  can  be 
found  by  moment  distribution.   The  deflections  at  each  cable 
attachment  of  the  frame  can  be  found  by  superposing  all  the 
deflections  at  that  point,  that  is,  the  deflection  at  any 
cable  attachment  of  the  structure  is  the  summation  of  deflec- 
tions on  that  position  due  to  external  loads  plus  the  unit 
deflection  on  that  position  due  to  unit  load  on  that  cable 
times  the  cable  stress  plus  the  unit  deflection  on  that  posi- 
tion due  to  unit  load  on  the  cable  at  other  position  times 
the  cable  stress.   Therefore,  for  each  redundant,  there  is 
one  linear  equation  of  superposed  deflections  to  be  written. 
For  many  redundants,  there  are  many  simultaneous  linear 
equations  to  be  solved.   The  results  of  those  simultaneous 
equations  are  the  cable  stresses. 

The  second  method  is  to  solve  the  cable  stiffened  struc- 
ture by  means  of  haunch  member  theory,  because  in  this  long 
span  cable  stiffened  rigid  frame,  those  cables  can  be  regarded 
as  part  of  the  frame  reinforcement.   By  the  method  of  trans- 
formed area,  those  cables  can  be  transfered  to  be  fictitious 
concrete  sectional  area.   Therefore,  the  uniform  sectional 
rigid  frame  can,  be  treated  as  a  frame  with  variable  sections 
without  increasing  concrete  quantities.   In  this  study,  the 
horizontal  girder  of  the  frame  is  first  divided  into  many 
equal  interval  sections  and  then  the  center  of  gravity  is 
located  for  each  section  of  the  concrete  with  transformed 


cable  sections.  The  mement  of  inertias  can  be  found  for  each 
section  in  respect  to  its  central  axis.   So  a  variable  moment 
of  inertia  can  be  found  due  to  the  variable  position  and 
quantities  of  cables.   The  next  step  is  to  solve  the  fixed 
end  moments  for  the  fictitious  variable  sectional  rigid  frame. 
The  theory  and  numerical  example  has  been  shown  in  this  re- 
port.  In  solving  those  Laborious  calculations  for  each  sec- 
tion and  for  the  fixed  end  moment  influence  lines,  a  FORGO 
program  has  been  provided  for  solving  this  type  of  problem. 

The  purpose  of  this  report  is  to  show  how  a  structure  can 
be  stiffened  by  cables  and  what  is  their  advantage.   Larger 
space  and  lighter  weight  of  construction  is  needed  in  modern 
architectural  design.   Cable  stiffening  will  enable  this  kind 
of  structure  to  be  built.   The  cable  stiffened  rigid  frame 
can  be  designed  for  many  purposes,  such  as  the  construction 
of  a  large  spaced  field  house,  auditorium,  factory  and  hangar. 
For  the  limitation  of  time  and  pages,  the  design  of  cable 
stiffened  structure  will  not  be  included  in  this  report. 
The  author  may  continue  further  study  after  finishing  this 
report . 


ELASTIC  PROPERTIES  OF  CABLES 

Cable  is  used  for  structural  purposes  in  a  variety  of 
ways,  such  as  diagonal  bracing  in  certain  types  of  frames, 
suspended  structural  system,  and  as  a  stiffening  for  some 
structures  for  the  purpose  of  enlarging  the  structural  space 
without  too  much  increasing  in  sectional  area  and  reinforce- 
ment . 

Since  a  cable  is  very  flexible  and  cannot  resist  bending 
or  compression,  the  characteristic  which  governs  its  effect 
upon  the  elastic  behavior  of  a  structure  into  which  it  is 
incorporated  is  the  product  of  its  cross  section  area  and  the 
value  of  Young's  modulus  in  tension,  AE,  which  will  be  referred 
to  as  the  extensional  rigidity. 

For  a  straight  bar,  the  extensional  rigidity  AE  is 
readily  determinable,  but  if  a  cable  sagging  under  its  own 
weight  is  subject  to  a  pull  tending  to  straighten  it,  the 
increase  in  distance  between  the  two  ends  depends  not  only 
upon  the  elastic  stretch  of  the  cable  but  also  upon  the  alter- 
ation in  its  configuration  due  to  the  change  in  the  dip-span 
ratio.   The  effective  strain  will  thus  be  greater  than  that 
calculated  from  the  elastic  extension  of  the  cable.   In  Fig.  1, 
a  cable  AB  of  length  L  is  hanged  at  point  A  and  B,  the  ver- 
tical distance  is  h  and  the  horizontal  distance  is  d,  and 
T  denotes  the  tension  at  each  end.   If  point  A  can  be  displaced 
horizontally  by  a  small  amount  Sd,    the  tension  at  any  point 


is  increased  by  £t.   Let  AE  be  the 

extensional  rigidity  for  the  material 
of  the  cable  determined  from  a 
straight  length,  and  7JAE,  the 
effective  extensional  rigidity 
of  the  equivalent  straight 
cable  required  to  replace 
the  actual  cable  AB.   Then 


oL  •=  £dcos  6  ■ 


L 


>[kl 


S*n 


^         AEd  S  a 


(1) 


where  T,,  denote  the  tension  in  the  equivalent  cable.   Let 
the  cable  be  a  constant  cross  section  and  total  weight  W.   It 
will  then  hang  in  a  column,  but  the  loss  of  accuracy  in  as- 
suming the  curve  to  be  parabolic  is  negligible.   Take  the 
vertex  of  the  parabola  as  origin  and  let  the  x-coordinates 
of  A  and  B  be  Xljd  respectively,   The  equation  for  the 
parabola  is  then, 

X2  =  2ay 

The  constant  horizontal  tension  is 


where  w  =  W/d 


also  T2  =  w2(a2 
and  X  =  atan  6 


(2) 
(3) 
(4) 


Where  a   is  given  by 


a        -sin  9  +7A(Ta/W)2  -   cos2  6 
d  2   sec   G 


sin  6  +/7UTB//W)      -    cos^e 
2   sec   9 


(5) 


Since  cos  o   is  small  in  comparison  with  any  reasonable  value 
of  T/W,  to  a  close  approximation  from  (5) 


TA  "  TB  =  Wsin  9 


(T.  t  TR)  -  HsecS  -  TV 


XA  T  XB 


(6) 


From  (3)  and  (4)  T  is  the  actual  tension  at  the  point 
x=X,  the  length  of  the  cable  is 


l   f (X+Jd)  T,    (X-sd)   T 


X*Jk   ^ 


X.-B+  alog( f- 


(7) 


A  small  change  <£d  in  d  if  h  remains  contant,  gives, 

From    (3)      %  .  dlM.xJx+  a^a  (8) 

WW 

From    (4)      h  <£  a  =  d£x  +  X<$d  (9) 

From    (7)     £s  =  A  £  X  +  B  S  d  +  C  /a  (10) 


Where  A  =  -A_ — £ 


TA  +  Ti 
2H 


10 

C   -  i(S-AX-Bd) 
a 

(11) 

Also     JS   =  -|^T 
AE        M 

(12) 

- 

Eliminating    £x  and^a  from    IB)    (9) 

(10)   and   substituting 

for     &S   from    (12) 

w     (^T(Ah+Cd)   "  l|(hX+ad))    =  ,Sd(AaX-B(hX+ad)    -CX2j       (13) 

The  two   terms   on  the   left   hand   side 

are   comparable   in  size, 

Ah+Cd   and   l/AE  both   small. 

By  using  the  approximation    (6) 

putting  S=L 

i(hX  +  ad)   -   L2Jf 

AaX   -   B(hX  +  ad)    -   CX2  =  -d^M 

(14) 

From    (11)   AhtCd  -  S   -    (TAI^B)    -   d  = 
da              2H            a 

/X+d 

i(       V(a2+x2)dx" 

f7a2+(X+id)2  +ya2+(X-id)2) 

Put  x'    =  x-   atan9,    than  by    (4) 

Ad 

Ah+Cd  „   secBf 1         /       x'(2atan9 
d              a     i      ,/            a2sec26 

+x'   dx' 

1  ,  1.  ,  d  (d-t^atane)     ,     /,  ,  d(d-4atane."l 
4a .'•sec^Q                     4a^sec^B     ) 

11 


On  expanding  the  right  hand  side  and  integrating  the  expression 
reduces  to  the  converging  series: 

Ah+Cd     l/dcoso,3    (l-5sin28 ) /deosg ^ 

a    "  12l   a  '      '  155    [T~'  (15) 

Retaining  the  first  term  of  this  series  and  substituting  from 
this  and  (14)  in  (13)  that, 

oTr.T  m  y 

"^  "  1     ^cos3fl+wL2  (16) 

n  — cos  9+ae 

Kence   from    (l) 

n= — srj (17) 

1+  ^Afcos2e 

'        12Tf3 
Also  from    (16)  , 

12T3 

and  the  total  displacement  of  A  as  the  mean  tension  in  the 
cable  increases  from  TQ  to  T,T  is 

A=  Lsecfi  fT  _T  i  ,  W_d   ,1    i  .  ,   . 

AS   UK  V  +  ZL        [W  ~  W]  (19) 

0    M 

Since  ^  depende  upon  the  force  in  the  member,  it  has  a 
different  value  for  every  loading  condition  of  the  structure. 
This  variation  ought  theoretically  to  be  taken  into  account 


12 


in  calculation,  but  this  would  be  impracticable.   If  the  load 
is  applied  in  increments,  }[   being  assumed  constant  during 
each  stage,  any  required  degree  of  accuracy  can  be  obtained. 
To  illustrate  the  importance  of  accurate  control  of  the 
initial  tension  of  cable  stiffened  structures,  the  case  of  a 
beam  hanging  at  four  equal  points  by  flexible  cables  will  be 
considered  in  some  detail.   The  beam  shown  in  Fig.  2  is  a  beam 
of  flexural  rigidity  EI  pinned  at  0,  and  hanged  at  A,  B,  C,  D 
by  cables  of  equal  equivalent  extensional  rigidity  AE,  with  a 
load  w  of  uniform  intensity  along  the  beam. 


Pig.  2 


There  are  three  redundancies  which  will  be  taken  to  be  the 
cables  FD ,  FC ,  and  FB.   Let  the  loads  in  these  cables  be 
RSSecQ,,  R2Sec92,  RlSecB^  respectively,  where  Rl,  R2,  R3  are 
the  vertical  component  loads  acting  perpendicular  to  the  beam. 
Then  by  taking  moment  about  0,  the  corresponding  component  of 
FA  is  SwL  -  (2R1  +  3R2  +  4R3 ) •   The  total  strain  energy  in 
the  structure  consists  of  that  due  to  the  bending  of  the  beam 
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and  that  due  to  the 

extension  of  the  cables  and  will 

be , 

2>U  -  2>U  „  3U 
2>R1   2>R2   3R3 

-  0 

The  contributions  to 

these  terms  from  the  bending  of 

the  beam 

are  deduced  from  the 

bending  moment  equations, 

2L.  =  il(  2,,,  . 

2R1   ET  3K1  + 

|R2  +  ?R3  -  ^wL) 

0R2    El  v  2ui 

4R2  +  2fR3  -  igwL  ) 

2>U   .  \? ,    7R,  + 

2R3  "  El(  3R1 

^■R2  +  12R3  -  ^fwL  ) 

The  components  from 

the  cables  are  found  as  follows, 

taking 

■2U/3R1  as  example, 

it  is  seen  that  the  unknown  Rl  1 

occurs 

in  the  cable  FA  and 

FB.   The  load  in  FA  is  P  =  (8wL  • 

-  2R1 

-  3R2  -  4R3)Sec6Q,  and 

m   i      .  pi 

13SI  J  FA   ^AE 

~^r   ,  where  I  denotes  the  length 

of  FA, 

that  is  hSec6  ,  also 

IL.  ■  -2Sec  6n 
3R1         ° 

v,^„,=  f  3U   i    _.  2h(2Rl+3R2+i-R3-8wL)Sec  60 
nenCS  I  SEE  J  FA               AE 

also  for  FB 

14 


C  au    "I  hniscc  a] 


>^AE 


[?et]      v.i        =-  -Ti7  2(2RH3R2+4R3-SwL)Soc3e+RlSec3e 
L  3fil  J  cables       rjAi, 


and  similarly   for   the  values   of  gSL  and  f&r   , 

2K2      ^'t-> 


Hence  the  total  value  of  M_  is 

?R1 


Sr  -  ^(fHi-^a^-ilwD  ♦ 


F|UURl+6R2+8R3-l6wL)Sec6+RlSec3e  J  -  0 

The  other  equations  are  obtained  similarly  and  finally, 

SRl+l8R2+28R3-39wL+12k  [(4Rl+6R2+8R3-l6wL)Sec3eo+RlSec3e1)  -  0 
l8Rl+4SR2+80R3-103wL+12k  [(6Rl+9R2+12R3-24wL)Sec3eo-fR2Sec3e2)  =0 
28Rl+80R2+144R3-174wL-H2k[(8Rl+12R2+l6R3-32wL)Sec3eo+R3Sec3e3)=0 

where  k  -  -&-  SI  • 
XJAE  L3 

Two  limiting  cases  provide  a  check  on  these  equations.   In 
the  first  case  if  fjAE  is  infinitely  large,  the  beam  cannot 
deflect  at  the  points  of  attachment  of  the  cables  and  becomes 
a  beam  on  a  series  of  fixed  supports,  and  can  then  be  analyzed 
by  the  theorem  of  three  moments, 
k  =  0,  and  the  equation  reduce  to 
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8R1  +18R2  +28R3  -  39wL  -  0 
16R1  -t48R2  +8QR3  -103wL  -  0 
28R1  i80R2  +14AR3-174wL  =  0 

The  solution  of  these  equations  is 

PA  =  R2  -  8/7  wL 

R3  -  13/14  wL, 

where  Pq  and  ?A  are  the  normal  forces  at  0  and  A.   These  forces 
are  the  same  as  those  obtained  by  applying  the  theorem  of 
three  moments. 

The  second  check  is  offorded  by  the  case  when  EI  is 
infinite,  that  is  the  beam  is  rigid,   k  is  then  ct>  and  the 
equations  become 


(4R1  +6R2  +8R3  -l6wL)Sec36Q  +  RlSec3e]_  =  0 
(6R1  +9R2  +12R3-24wL)Sec;3e0  4  R2Sec3G2  -  0 
(8R1  +12R2+l6R3-32wL)Sec3Gn  +  R3Sec3e,  =  0 


The  solution  of  these  equations  is 

3 
SwLCos  6 
PA  ** 


Cos36qh  4Cos3b1+9Cos3e2+l6Cos3e3 


R  _  l6wLCos3ei 


Cos3d0+4Cos3eii9Cos3e  +l6Cos39, 
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R2 


JZiwLCosJt)~ 


Cos3e0+i|.Cos3e14  9Cos  8~-\  l6Cos  9o 


32wLCos36 


R3  " 


1 


3       3       3        3 
Cos^6qh  4Cos  e^Ccis  62h  l6Cos  63 


These  results  can  also  be  obtainable  from  a  direct  consider- 
ation of  the  displacements  of  the  rigid  beam. 

The  bending  moment  diagrams  resulting  from  these  two 
assumptions  are  shown  in  Fig.  3.   It  is  evident  that  the 
differences  are  considerable  but  neither  is  correct,  because 
neither  Y[AK   nor  EI  can  actually  be  infinite,  and  k  will  have 
a  finite  value  depending  upon  the  initial  tensions  in  the 
cables.   The  bending  moment  curves  are  shown  for  values  of 
0.01,  0.05,  0.1  and  1.0.   The  critical  point  for  the  design 
occurs  at  A,  the  bending  moments  there  are  given  in  the 
following  table  for  different  assumptions: 
k       0    0.01  0.05   0.1    1.0 

k=0  k=0.lO 


%  0.107      0.113   0.171     0.215      0.382      0.499 


0  A  B  C  D 

:c=0.01 

u  A  Ji  C  D 

k=0.05 

0  A  3  G  I) 


BCD 
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From  Fig.  3  the  very  big  differences  in  M.  emphasize  the  im- 
portance not  only  of  proper  control  of  initial  tensions  but 
also  of  accurate  methods  of  calculation. 

The  foregoing  analysis  has  been  simplified  by  the  re- 
striction that  all  cables  shall  have  the  same  extensional 
rigidity,  it  has  provided  no  means  whereby  the  beam  can  be 
maintained  in  a  straight  line.   It  serves  to  show  that  in  a 
condition  of  loading  in  which  the  cable  points  on  the  beam  do 
not  remain  collinear,  the  relative  elasticity  of  the  beam 
and  cables  becomes  an  important  factor.   It  is  assumed  that 
the  far  cable  ED  from  the  end  0  is  displaced  by  a  specified 
amount  and  that  the  cable  points  remain  collinear. 

From  this  assumption,  a  further  discussion  on  the  geom- 
etry of  the  vertical  deflection  of  the  beam  and  displacement 
of  cable  attachment  is  required.   In  Fig.  I,  OB  is  a  column, 
and  a  cable  BA  hangs  the  beam  OA  at  point  A,  For  the  sake  of 


Pig.  4 
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simplicity,  assume  the  beam  is  hinged  to  the  column  at  point 
0,  that  is,  there  is  no  fixed  end  moment  at  the  connection. 
A  small  deflection  is  assumed  on  this  structure,  that  is,  the 
angle  between  the  cable  and  column  does  not  change  as  the 
deformation  take  place,  and  the  horizontal  deformation  of 
point  A  is  negligible.   Now  the  same  geometrical  relations  can 
be  deduced  as  follows:   Let  ok   be  the  vertical  displacement 
at  point  A,  assume  downward  displacement  is  positive,  and  let 
Ax  be  the  horizontal  displacement  of  point  B   of  the  column, 
a   «  the  angle  between  the  cable  and  the  column,  L  =  the  orig- 
inal length  of  the  cable,  AL  »  elongation  of  the  cable,  As  = 
the  cross  section  area  of  the  cable,  Es  -  the  elasticity 
modulus  of  the  cable,  T  =  the  tension  in  the  cable, 

AL  -  B1A1  -  3A     since  B1A1  =  B]A2  +  A^ 

BA  =  BB,  +  B,A 

and  B,A  =  B-jA-, 

The  elongation  of  the  cable  BA  becomes, 

AL  -  B1A2  +  A2A1  -  BB3  -  RjAg  =  A^  -  BB 

"  Skcosd  -     xsin6  (20) 

AL  -  TL/EAs 

.'.   AcosS  =  -££      £k   =  ?-i_  (21) 

EAS  EAEcos9 
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The  vertical  cable  tension 

T,  -  JjMs    oos2e 

Letting  £A  ■  unity,  and  the  vertical  cable  tension  shall  be 

T>  =-  ^-s  cos2e                             (22) 

' 
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STRUCTURAL  ANALYSIS  BY  THE  PRINCIPLE  OF  SUPERPOSITION 

The  principle  of  superposition  applied  for  solving  high 
degree  indeterminate  structures  has  been  considered  the  most 
efficient  method  of  structural  analysis  by  using  a  Computer. 
In  applying  this  method  for  stress  analysis,  it  is  necessary 
to  reduce  the  structure  to  be  a  base  structure  form,  which  is 
a  statically  determinate  or  the  minimum  degree  statically 
indeterminate  structure,  and  put  all  the  redundants  to  be 
solved  by  linear  simultaneous  equations.   Those  equations  can 
be  solved  by  matrix  operation. 

The  principle  of  superposition  is  also  suitable  for 
analysis  of  cable  stiffened  structures.   Consider  the  cable 
stiffened  structure  as  shown  in  Fig.  5,   It  is  a  single  bent 
rigid  frame  with  six  cables  hanging  on  the  top  of  the  girder. 
It  is  a  statically  indeterminate  structure  in  nine  degrees. 
Ey  the  principle  of  superposition,  the  total  moment  distribu- 
tion can  be  considered  as  composed  of  two  parts, 

(1)  The  moment  distribution  mg  due  to  the  applied  loads 

acting  only  on  the  base  structure  without  cable  stiffening, 

that  is  the  cable  stresses  Xt=X„  .  .  .  .  =X  =  0.   It  is 

J-  2  n 

called  the  particular  solution  of  the  problem  because  it 
satisfies  the  condition  of  equilibrium  but  does  not  satisfy 
the  boundary  condition  of  the  problem,  since  there  are  six  X 
unknowns  to  be  solved. 

(2)  The  moment  distribution  m  X  due  to  the  action  of 

n  n 

the  arbitrary  constant  Xn  alone  on  the  structure  with  all  other 
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loads  equals  to  zero.   The  number  of  arbitrary  constants  it 
has  should  be  the  same  as  the  number  of  moment  distributions 
due  to  these  constants. 

Thus  the  total  moment  distributions  should  be 


M  -  m  +  m1X1  +  rruX  H  m,X 


m  X 
n  n 


(23) 


The  total  strain  energy  due  to  flexural  bending  as  discussed 
in  general  strength  of  materials  will  be 


U-  (j£ds 
2EI 


(24) 


It  is  also  necessary  to  consider  the  effect  of  stress  resul- 
tants other  than  moment  in  contributing  to  the  strain  energy, 
and  the  total  energy  will  be  the  summation  of  these  energies, 


total    bending   shear   direct  force 


In  equation  (24)   s  is  in  term  of  the  coordinate  and  is  as- 
sumed along  the  axes  of  the  members  of  a  skeletal  frame. 


A   X5  "X* 


^A 
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Applying  the  theorem  of  least  work,   n  equations  can  be 
obtained  for  the  determination  of  nX  values. 


wi  ~    \  \  {m)as  =  ° 


2U      (     ~d     ,   M2,, 


f*J  ~     J  X^  (iM)ds  =  °  (25) 


2U  _   /  3  ,  M2,  ,  _  n 
P^  "   J  X~n    (2EI)dS  -  ° 


It  is  known  that 

U  =  |  —ds   ■    -i-(mn+m,X,+m0X„+m,X, .  .  .m  X  )ds 
1  2EI     I  2EI  0  1  1  2  2  3  3      n  n' 


So  that  equation  (25)  becomes: 


|E-  =       §l(m+m  Y.,+m,X0+m.X,.  .    .m  X    )ds  =  0 
^X.  Ei      0      11      22      33  nn 


%T2  '       1f(VmlXl+m2X2+m3X3-    •    -mnXn)ds  "  ° 
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gU    -    I  ^(m  +m.X.+m0X0+m,X 


2>X,  KT    0      11      2   2      3   3 


,+m  X    )ds  -  0 
n  n 


2Xn         j   hi      0      11      22      33 


.    ,+m  X    )ds  -  0 

n  n 


It  is  convenient  to  expand  these  equations  as  follows: 


Xl  )  EIds  "•  X2    -Tru»   TA3  |   ""EI 
Js  Js  '   Js 


mn  in,  nil  m,  I     m-i  mn 

-?ds  +X,    \     -Mis.    .    .+X„    I     4^ds 


n  EI 

/s 


EI 


-ds  -  0 


m,mn  mi  /     mm 

.  )s  -TT^  +  X2  Js  lfds  +  x3   I    "lTldS-    • 


m9mr 


,  )s  -frds  +  Js  Tds "  ° 


m0mr 


,    j^ds   +  X2   • 


Sds- 

Js 


m,m  /     m,mn 

;  )s  ~irds  +  L  ~irds  '  ° 


mm.,  /     m  m„  /mm, 

s  4rds  +  x2  I  "ir*8  +  x3  (  -ii1^- 


m  (mm,. 


>      (26) 
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From  Castigliano's  theorem,  it  is  known  that 

f^J  =  $    .   So  the  deflection  <£,  at  the  position  and  the 
direction  of  X.  should  be: 

2 
m-t  I   mn  ntn        I   mi  m^ 


K-h-h\  »■♦*     ^'M  ¥*•• 


s 


n  L  KI       |_  EI 


ds  +    -i^ds 


(27) 


For  simplicity  equation  (27)  can  be  defined  as: 

L2 

*1  =  Xl  JETds  =  Xlfll    due  t0  Xl  =  1' 

where  f  .  =  the  deflection  of  the  base  structure  at  the 
position  and  in  the  direction  of  X-,  due  to  unit  X-,  acting 
alone. 

K   =  X2  ^^   -  Vl2    due  toX2  =  l, 

where  f -,  ?  =  the  deflection  of  the  base  structure  at  the 
position  and  in  the  direction  of  X  due  to  a  unit  X2  acting 
alone.   And  so  do  the  deflections  due  to  other  unit  X  acting 
alone  as, 

r  (  m.  m 

dl   =  X3J^-ds  =  Vl3     due  to  X3  =  1, 

r  [mm 

°1  "  Xn  Hr1  ds  -  V in    duetoXn  =  l, 
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C            (mlra0  . 

*1    J  KI   dE  °  ul 

due 

to  applied  loads, where  u-^  =  the 

deflection 

of  the 

base 

structure  at  the  position  and  in  the  direction 

of  Xx   ( 

iue  to 

the 

applied  loads. 

Now  with  the  above  definition, 

equation  (26)  can 

be  re- 

wri' 

;t en  in  the  following  form,  for 

the  first  one: 

Xlfll+X2f12+Vl3-  •  -+Xnfln= 

"ul 

(28a) 

Sim: 

Llarly,  the  second  equation,  as 

well  as  the 

rest  o: 

f  the 

equations  in  relation  to  each  cable 

as : 

:  attachment 

can  be 

written 

Xlf21+X2f22+)3f23+"  -  -+Xnf2n= 

:  "u2 

Xlf31+X2f32+X3f33+•  -  -+Xnf3n= 

■  -U3 

X,f  ,+X„£  „+X,f  ,.  .  ,+X  f  = 

1  nl  2  n2  3  n3       n  nn 

-u 
n 

(28b) 

Where  fnl,  fn2,  f..  .  .f   have  meaning  similar  to  f. 

LI'  f12' 

f13, 

.  .  .£,   and  u  also  have  the  same  meaning  as  u,  but  those 

In      n 

are 

related  to  Xn.   It  is  to  be  noted  that  f,2 

=  f21> 

f31  = 

f13 

•  "  -fnl  =  fln  because, 

t    m,m? 
f12  -  f^  -     ^"ds,   which 

Lprocal  theorem  due  to  the  physi 

is  from  the 

Maxwell's 

rec: 

cal  meaning 

of  f12 

and  f  „,  . 
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In  summation  of  the  above  statements,  a  general  solution 

of  a  structure  with  n  arbitary  constants  will  lead  to  a  set 

of  n  simultaneous  equations  of  the  form  of  equation  (28). 

The  equation  to  be  solved  will  then  be: 


Xlfll+X2f12+):'3f13' 


•VU-  •  -Via  =  "Ul 


1  21  2  22  3  23      J  2j      n  2n     2 


x1f31+x2f32+x3f33.  .  .Xjf3j-.  .  .xnf3n  -  -u3 


XlfirX2fi2+X3fi3- 


.X.f X  f.   =  -u 

J  ij      n  in     l 


>      (29) 


X.f  -,+X„f  „+X,f  ,  .  .  .X.f  ..  .  .X  f   =  -u 
1  nl  4  n2  3  n3      j  nj      n  nn     n      J 


It  should  be  noted  that  f . .  ,  f,..  £-,■,.    .    .f   are  the  de- 
ll'  22'   33      nn 

flections  of  the  base  structure  at  the  position  of  cable 
attachment  due  to  a  unit  load  acting  at  this  point,  and  the 
actual  deflection  at  these  points  should  include  the  elonga- 
tion of  the  cable  at  that  point,  that  is 


1, 


"1    f 


2  I    9 

n,      L„         {mf     L_ 

AA  c  -+■  fc 


hi  A^E'    22   JirS  +  A^E'  '••   nn  JlTds  IT! 


Where  f i  .  is  called  the  influence  coefficients,  which  is  a 
function  of  X.   Equation  (29)  can  be  written  in  matrix  form 
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f     f     f 
11   rl2   *13 


f21  f.,.,  f. 


"22    23 


f31   f32   f33 


f     f     f 
il    12   ri3 


f     f     f 
I  nl    n2   In3 


fU 


f2J  • 


f 
3j 


iJ 


nj 


In 

'  V 

ul 

2n 

X2 

u 
2 

3n 

X 

3 

^~ 

u 
3 

in 

X. 

J 

u. 

i 

nn 

X 
n 

u 
n 

Usually,  it  can"  be  written  in  a  shorthand  notation  as 

[F)(x)=-(u)  (30) 

Here,  (fJ  is  called  Flexibility  matrix. 

The  solution  of  the  equations  can  be  obtained  by  matrix 
operation.   Usually,  this  matrix  operation  can  be  transferred 
into  mechanical  language  and  solved  by  digital  computer.   If 
it  is  solved  by  hand,  the  following  operation  is  suggested. 

First,  premultiply  each  side  of  equation  (30)  by  F  _1, 
that  i  s 

(F)-1(F](X)=  -(F^lu) 

Since  [F]  _1  [F]  =  [I]  is  an  Idendity  matrix,  therefore,  the 
solution  can  be  written  in  the  form  as 

[X]  =  -[F]"1  [U]  (31) 

where  ( F )    is  the  inversion  of  matrix  [i\.      It  is  obvious 
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that  any  kind  of  high  degree  statically  indeterminate 
structures  can  he  analysed  by  matrix  method,  and  it  is  time 
saving  and  convenient  to  solve  this  problem  by  computer. 

In  addition,  the  author  wants  to  introduce  some  matrix 
operating  concerning  this  kind  of  problem.   From  equation  (30) 
since  it  is  known  that  (f)  is  a  symmetrical  matrix,  then  there 
can  be  found  a  lower  triangle  matrix  [L],  such  that 

where  (l)'  is  the  transpose  of  matrix  (L),  which  is  an  upper 
triangular  matrix.  Now  the  operation  of  expression  can  be 
carried  out  in  fact  as 


Lll 

0 

0   .  . 

.  0 

L21 

L22 

0   .  . 

.  0 

L31 

L32 

L33  •  • 

.   0 

Lnl 


Ln2 


Ln3 


Lnn 


Lll  L21  L31 
0  L22  L23 
0    0   L33 


0 


0 


0 


Lnl 
Ln2 

Ln3 

Lnn 


Lll    L12    LI 3  .  .  .  Lin] 

L21   L22   L23  .  .  .  L2n 

L31   L32   L33  .  .  .  L3n 

. Lnl   Ln2   Ln3  .  .  .  Lnn , 

So  from  the  above  expression,  equation  (30)  can  be  written: 

[LJU)'(X)-  -£u]  (32) 

Similarly,  [F]"1  =  ((L)  [L  ]' )_1  =  (l)-1(l)'-1 
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Postmultply  each  side  by  (l]  ,  and  the  above  equation  becomes 

Therefore,  in  matrix  inversion,  it  will  be  easier  to  find  [  L  ) 
and  (L')~  than  to  find  (  f]_1  directly. 

Two  numerical  examples  will  be  given  as  follows  to 
demonstrate  how  to  apply  the  principle  of  superposition  to 
solve  the  cable  stiffened  structural  problems. 
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NUMERICAL  EXAMPLE  I 

Analyze  and  design  the  stresses  of  a  cable  stiffened  rigid 
frame  as  shown  in  Fig.  6.  The  frame  has  a  span  of  140  ft. 
center  to  center  of  columns  and  the  lower  column  is  50  ft.  in 
height  and  the  upper  column  is  20  ft.   The  section  of  the  beam 
has  5  ft.  depth  with  2  ft.  width  and  the  section  of  columns 
is  12  ft.  depth  by  2  ft  width  which  is  fixed  supported  at  each 
end.   The  frame  is  stiffened  by  six  #lSs  high  tension  steel 
bars  as  cables  hanging  on  the  top  of  the  beam  and  with  equal 
distance  of  20  ft.  apart.   Space  between  each  frame  is  10  ft. 
with  4"  roof  slab  on  the  top  and  12"  wall  on  each  side  of  the 
frame .   Assume 

fg  =  35,000  psi  for  tension  cable  bars 

f  =  60,000  psi,   fs  =  24,000  psi  for  frame  reinforcement 


"12"  wall 


1  =  2oo. 
c 


Ib=20.833 


Pig.  6  Frame  Layout 


-liiO-'- 


V  slab 


50  > 
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fc  -  3 ,000  psi 

.\.  =■  29,000,000  psi  1SC  -  2,500,000  psi 

Cable  length  LI  -  23.235'   L2  -  44.745'   L3  -  63.253' 

Loading  condition 

Dead  load   beam  5x2x0.15        =  1.5  k/ft 
slab  (4/12)3x0.15     =0.4  k/ft 
Live  load   snow  35/v/sq.ft  x  10'   =  0.35/ft 
The  wind  stresses  will  be  shown  in  the  next  example 

total  load  =2.25  k/ft 

Check  of  cable  size 

Use  of  #l8a  as  tension  cable  (2.257")   As=4.0  sq.  in  = 

0.0277S  sq.  ft. 

AS  =   I  4-x  sec  e    25  x  5x2-25xl402xl.4142  ^ 
Nxf^      "        3x35000 

3.62  sq.  in  <  4.0    0  K 

Where  6  is  the  angle  between  the  cable  and  the  beam,  S~ax  = 
45°  for  cable  1,  h  is  the  length  of  upper  column,  K  is  the 
number  of  cables.   This  formula  is  derived  on  the  assumption 
that  the  bending  moment  of  the  beam  is  taken  by  those  cables. 
So  the  assumed  As  is  ok. 

Moment  for  the  base  structure: 

The  base  structure  is  shown  in  Fig.  7,  which  is  the  orig- 
inal frame  uniformly  loaded  of  2.25  k/ft  on  the  top  of  the 


2gie/ft 


tt,  4-  v  i  nrf  uihh  wan  tf 
1-2  v.  1.33  i't."' 


K=0.1l<9 
S'=0.5S=0.07ii5 


JM£L1 


(a) 
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beam  without  cable  stiffening.   The  fixed  end  moment  and  the 
moment  diagram  can  be  calculated  by  the  method  of  moment 
distribution.   (NOTE:   All  moments  are  in  kips-ft.  in  this 
report . ) 

Moment  for  the  base  structure  due  to  unit  tension  on  cable: 

The  base  structure  is  loaded  by  cable  1  with  a  tension 
force  X  =  1  on  the  cable  and  the  fixed  end  moment  is  found 
as  shown  in  Fig.  8a. 
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Korient   of  the  base   structure   due   to  unit   tension  on  cable   2: 
(see   Fig.    6b) 

..0.69395  0.69395 

0.41l696y. 


X=1      8.44696   «-H4p8x_1 


JC=5.76 


K'=  0.0745 
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Itowent   of  the  base   structure   due   to  unit  load  on  cable   3: 
(see  Fig.    6c) 

_0. 94857  0.94657      _ 
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(a)  Koruent  diagram  of  base  structure  due  to  external  loads 
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(b)   rloraent   diagram  of  base   structure   due  to  unit  load 
on  cable   1  . 
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(a)  Moment  diagram  of  base  structure  due  to  unit  load 
on  cable  2.      V 
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(b )  Moment  diagrara  of  base  structure  due  to  unit  load 
on  cable  3- 
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Deflections  of  the  base  structure: 

The  producing  and  integrating  of  the  moment  areas  have 
been  done  by  Computer,  Programs  and  data  output  for  this 
problem  are  shown  in  Appendix  B.  So  the  deflections  due  to 
external  loads  and  unit  loads  from  the  Computer  output  shall 
be: 


mQm1 
1   1  ~I2~  ds  =  -  18070. 64300/E 


mnm„ 
U2=  I  TT^  ds  =  -31524. 613/E 


-§jp.  ds  -  -31746. 20100/E 


/  m-  ^  L 

fll=  (  -gj  ds  +  — I  *  +55. 089719/E+1018. 1785=1073. 268219/E 

s  1 

m-imo 

f12&  i  -fr   ds  =  72.755847/E 


/  m.m, 
f13"  \     Til  ds  =  +66.19547/E 
s 


i  da  +  ±1  - 


22°    EI  ds  "*  A~E   125. 375590/E+1610. 6911  =  1736.066690/E 


f23=    ~fi  ds  "  123-3H750/E 
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-|j±  ds  -  72.755847/E 


J    He 


-£  ds  +  -A  "  134.565680/E+2276.9258/E=2411.49180/E 

jj  ai      ^"3 


m3m2 


ds  -  123.3H750/E 


Sh*    \    ~ET   ds  =  66-1955470/E 


Those  influence  coefficients  can  be  put  into  a  matrix  equation 
as: 

1073.268219   72.755847    66.195547' 
72.755847  1736.066690   123.3H750 
I  66.195547  123.314750  2411.491480  J 


rV 

'  18070.643 

X2 

31524.613 

U3 

,  31746.201 

The  inversion  of  the  (F)  matrix  and  the  value  for  (X)  matrix 
have  been  done  by  the  Computer  program  3-  So  the  cable 

stresses  are  found: 


[X] 


V 

X2 

= 

K 

IF]"1 

0.00093574      .  -0.000037527  -0.000023767 

-0.000037527        0.00057962  -0.000023767 

-0.000023767     -0.000028609  0. 00041680 


-luj 

-u. 


l-3, 


14.971933 
16.839623 
11.900315 


■S&.&4557  \ 


1 


V 

s 
I 


'0  ^ 


\   a    a 
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4 

I 

0         1   0.98?  j 

O.013 

I 

-732.91il;.u       0 

0       +4.96012 

+727.86569 

+0.06559 

f 

-732.911^+I|.9('012 

+727.93120 

+2.1J.9006 

(a)  Loads  on  frame  d-as  to  cable  stresses  and  moment 
distribution. 
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(b)  Komont  diagram  of  the  frame  .*»  t0  oable  stresses. 


Pig.  10 
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Fig.  11  Moment  Diagram  of  the  Swaua  due  to  Sxtornal  Load3 
and  Cable  stresses 


Fig.  11  shows  the  actual  moment  diagram  of  the  cable 
stiffened  rigid  frame  due  to  external  loads.   The  moment  at 
the  fixed  end  of  the  beam  is  2399.29372  kips-ft.  in  compara- 
tive to  the  fixed  end  moment  3627.2250  kips-ft.  of  the  beam 
without  cable  stiffening.   So,  the  moment  of  the  beam  has 
been  reduced  in  an  amount  of  727-93128  kips-ft.  due  to  the 
counter-balance  of  the  cable  stresses  with  the  external  loads 
on  the  beam. 

If  the  moment  of  2899.29372  kips-ft.  is  still  larger  for 
the  beam  to  take,  a  further  reduction  of  the  moment  can  be 
made  by  prestressing  the  cables.   The  amount  of  moment  reduc- 
tion desired  can  be  interpreted  as  the  amount  of  prestressing 
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of  each  cable  needed.   Assume  that  the  fixed  end  moment  of 
the  beam  need  to  be  further  reduced  in  a  amount  of  M  kips-ft., 
which  can  be  done  by  prestressing  those  cables.   The  cable 
prestressing  can  be  done  by  the  following  ways:   Figure  12 
shows  a  construction  detail  of  the  cable  post-tension  method. 
It  is  done  by  putting  a  small  jack  between  the  cables  and 

the  top  of  the 
concrete  filled  after  prestressing    columns  on  each 

side  of  the  frame. 
After  the  construc- 
tion of  the  frame 
is  finished  and 
has  reached  28 
days  concrete 
strength,  then  the 
cable  should  be 
jacked  up  in  a 
distance  Ah.   The  amount  of  Ah  can  be  determined  by  means  of 
interpreting  which  moments  should  be  reduced  by  cable  pre- 
stressing. 

It  is  known  from  equation  (21)  that  the  vertical  elonga- 
tion of  the  cable  shall  be 


Pig.  12  Cable  post-tensioning  Layout 


Ah  = 


TL 


EA  co; 
s 


where  T  is  the  tension  in  the  cable. 


and  the  vertical  component  of  the  cable  tension  shall  be 
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AhEA„ 


2 

cos  e 


where  6  is  the  angle  between  cable  and  column.   So  the  opposite 
fixed  moment  due  to  cable  prestressing  can  be  found  without 
difficulty.   Let  T'i,  T'0,  T'   be  the  vertical  components  of 
the  cable  tension,  and  the  corresponding  angles  between  column 
and  cable  1,  2,  3  be  6^  62,  9,  respectively.   Then  the  fixed 
end  moment  due  to  cable  prestressing  will  be: 


M  " 


T ' x ( 20x120^ 120x20^ )    I ' 2 ( 40x100^100x40^ ) 


140' 


HO' 


T'3(60x802+80x602) 
140  2 


•M  =  17.14290T'  +  28.57142T'  +  34-28571T' 


The  vertical  component  of  cable  prestressing  in  terms  of 
vertical  displacement  for  each  cable  will  be: 


T'   =  AhEAslcos2e,   T' 


AhEA 


1 


L. 


^As2C0S2ej      T,    =     AhEAs3C0S2e 


Vfeere   I*-   28.285'  L  =  44.745'        L,=  63.253' 


cos91  =  20/28.285  »  0.70708  cos291  =  0.49996 
cos92  =  20/44.745  =  0.44698  cos292  =  0.19979 
cosS,  =  20/63.253  -  0.31620   cos293  =  0.09998 
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Thus 

T'   -  0.000/, 91  E  h,   T'   =  0.000125  E  h, 

T'  =  0.0000439  E  h.  Substitute  T'  and  found  the 
moment  in  terms  of  E  h. 
Thus 

H  =  0.00S417  EAh  +  0.003571  E  Ah  4  0.0015219  EAh  = 
0.013510  EAh. 
Now  assume  that  the  fixed  end  moment  of  the  beam  will  be 
reduced  in  an  amount  of  1000  kips-ft.  so  the  prestressing 
should  be  done  by  jacking  up  the  cables  in  a  vertical  dis- 
tance Ah,  which  will  be  equal  to 

Ah  -   1000    „      1000x1000 ■>/-£  ,  ,~7ft 

0.013510E   0.013510x29,000,000x144    '*».«J.//o 

-  0.01776  ft. 

The  moment  diagram  due  to  cable  prestressing  and  the 
final  moment  and  shear  diagrams  are  shown  in  Fig.  13  and  Fig. 

14. 
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(a)  Loads  on  fpaine  duo  to  cnble  prestressing 
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(b)   Koment   diagram  du<j   to   cable  prestressing         w 
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(o)  Final  moment  diagram,  gf  the  rigid  frame 
Pig.  13 
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(a)  Shear  diagram  of  the  base  structure  due  to  external 

loads  without  cable  stiffening 


i/qr&4&£ 


1175/443 


(b)  Pinal  shear  diagram  of  the  frame  due  to  external  load 
and  cable  prestresoing. 
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NUMERICAL  EXAMPLE  2 

Analyze  the  cable  stiffened  rigid  frame  as  shown  in 
numerical  example  1,  loaded  by  lateral  wind  loads  of  20  lbs. 
per  sq.  ft.  on  the  left  side.   No  live  load  is  considered  on 
the  top  of  the  frame,  but  about  2  k/ft.  of  beam  and  slab 
weight  will  be  acting  on  the  top  of  the  frame. 

This  problem  can  be  analyzed  by  two  parts;  one  is  sym- 
metrical loading  condition  and  the  other  is  anti-symmetrical 
loading  condition.  The  summation  of  the  analyzed  results  of 
these  two  loading  conditions  will  give  the  actual  result  of 


ictual  condition   0.2" 


Ib=20.633^. 

I  =288 
c 


anti-symmetrical 


0.1 


i</rt 


K<  =0.2235' 

v=5.76 


+ 


-JD.1    k/ft. 

i 


Byrasaatrical 


0.1 


k/f 


K'=o.071|5  l 
K=5.76 


Fig.   15 


-0.1    k/ft. 


a 
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the  frame  stress  analysis  as  shown  in  Fig*  15.   The  lateral 

wind  stresses  on  the  frame  will  be  0.02  x  10  -  0.2  k/ft. 

Moment  for  the  base  structure  due  to  external  load  (symmetri- 

cal): 

5c/: 

O.v 

k   2ic/ft. 

■J   1  r   -t  *"1  v  i  r~i  *  t~±  i   I   I    I   1  1  1    1.  ] — r 

j 

0 

O.967 

0.013 

:G                 '1=20.633 

:  1=288         X=0.'i!i.9 
:K=5,76       K'=0.5K=0.07) 

I  k/ 

;0.1 

'-20   +20.8333 

0  +3232.377! 

-3266.666     L 
+^2. 1, 5583      I 

-20   +32I|1)..210( 

3  -32)|i|..2lu83 

+1611 .68675 

(a)                         -20.8333 

+1590.S551|5 

Moment  for  the  base  structure  due  to  external  loads 

( anti-symmetrical ) : 

k/4 

O.V 

-1 

k  °k/1     6 

0 

0.962 

0.036 

j 

^~- ^ " 

•  1=20.633    * 

1=286       K=o.ilt9     Z 

K=5.76       K'=1.5X     - 

=o.2?35   : 

-20   +20.8333 
0   -0.601 35 

-0.03165 

20   +20.03165 

-c. 031 65 

(b)                        -20.6333 

-0.U0678 

Pig.  16,  a  b                      _21 .23978 
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Moment  for  the  base  structure  due  to  unit  loads  (symmetrical): 

The  symmetrical  condition  of  unit  load  on  cables  had  been 
shown  in  example  1,  and  the  calculations  and  moment  diagrams 
can  be  used  in  this  example.   Moment  diagrams  refer  to  Fig.  9 
Co),    (c),  (d). 

Moment  for  the  base  structure  due  to  unit  loads  (anti- 
symmetrical)  : 

Anti-.sjmmetrioal  unit  load  on  cable  1  : 


0.70708       c. 70706  _~ 
A%f-  70705    o.707oCyx='' 


K-' =0.2235 


K=5.76 


+2.63763 


(a) 


Anti-synraetricul  unit  load  on  cable  2: 


0           0 

.962 

c. 

nd 

-1).!. 

.11*16 

0 

+6 

.65637 

\ 

0    +5., 

?7525 

+0 

20830 

-11. 

.11*16+5 

2752 

;+8 

86675 

.89395 
0.1*1* 

698 

0 

0.89395 
.1^695^, 

•> 

K=5 

76 

K'  = 

-^ 

.2235 

| 

0          0.962 

0.0  38 

-17.0790      0 
0    +11.931+3; 

-17. 8790+11.  931|0 

+5.U7323         L 

:+o. 1*711+2             "f 
,5+5.91(465 

+5.96716 


(b) 
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Anti-symmetrical  unit  load  on  cable  3: 


r 


y.<9U6<"! 


o. 94657. 


XH  °-31(M9Oh3l6X=1 

raJ L^ 


o 


0.962 


-18.971U0   0 
0   +16.76O6 


-18.9714+1 6. 7606 


0.038 


n .54869 

:+o.  66206 


+2.21075 


+8. 38033 


(c) 


Pig.  17 

The  moment  diagrams  both  for  external  load  and  unit  loads 
on  the  base  structure  are  shown  in  Fig.  18  and  Fig.  19. 

Deflections  of  the  base  structure: 


The  multiplying  and  integrating  of  the  moment  areas  have 
beer,  done  by  Computer,  and  the  program  and  data  output  for 
this  problem  is  shown  in  Appendix  B,  the  program  2.   The 
inversion  of  the  influence  coefficient  matrix  f.  .  and  the 
final  solution  of  the  cable  stresses  X-L,  X.,   X,  due  to  wind 
load  and  dead  load  of  the  frame  are  also  done  by  the  Computer; 
the  answer  is  shown  in  Appendix  B,  program  3.   and  the  output 
of  problem  2. 
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(m0  diagram) 
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(a)   I-jorncnt   diagram  of  base   structure  due   to   external 
loads    (symmetrical   loading   condition) 
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(b)  Moment  diagram  of  the  base  structure  due  to 

external  loads  (anti-symmetrical  loadind  condition) 
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(a)  KoJttent  diagram  of  base  structure  due  to  unit  load 
on  cable  1.  (anti-symmetrical  condition) 

V 


sa&rts  S9en&' 

(b)  ;  o.iont  diagram  of  base  structure -due  to  unit  load 
on  cable  2. (anti symmetrical  condition) 
Fig.  19 
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la)  Moment  diagram  due  to  unit  load  on  cable  3« 

Kg.  19  continue 
Deflections  of  the  base  structure:   (symmetrical  load  con- 
dition) 


HJT'ds  -  -11776. 1340/B 


u0m2 


=^ds  =  -23875. 7760/E 


3  \     ~Yp-ds  =  -27264. 7400/E 


lll 


L12 


\    ~E^ds  +  r^  =  35.60180+  1018.1785  -  1053.78030/E 

)s  1 


4f^ds  =  56.02362/E 


: 


t      m,m, 
LI3  J     -gf^-ds  -  51.87099/E 
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1       1       EI 


kis  -  56.02362/E 


22=    \    Elds  +  A~E       107-a3205+l6l0.6911  -  1718.52315/E 


V»3 


^ds  =  107.27473/E 
f31°  \   Eiids  =  51.87090/E 

f32=  (  ^|x2ds  "  107.27473/E 
/  mm       L, 

33~  \    ei  ds  +  ri=  129.05937+2276.9258  =  2405.98517/E 

/s  3 

Those  influence  coefficients  can  be  put  into  a  matrix  equation 


1053.78030    56.02362    51.87099 
56.02362  1718. 52315   107.27478 


-11776.13400 
-23875.77600 


.    51.87099     107.27478    2405.98517  J  [X, J      [-27264.74000, 

The  inversion  of  the  matrix  and  the  solution  for  the  cable 
stressed  have  been  done  by  the  Computer  (see  Appendix  b), 
answer  as  follows: 
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^3 
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^   esr4irr 

I 


symmetrical   condition 


(a)   Loads   on  the   frame   due   to   cable   stresses 


■  0 


o.9e? 


■574.835i(.o      o 

0 J-1  70 .11,2(119 


-5.71) . 633'J+-iYi  70.i(.20V9 


5ty6%4 


V, 

O 

6S 

vS 

t 

<Ti 

1 

8 


a 


A  ~\ 


+402. 17057 

■l-2.2iil|hli-     - 


+I|.0i|..l|l521 


V  +85. 2^010 

$  '   ft 


£5.81010 
(b )    .  .oaisnt  diagram  due   to   cable   stresses 
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V 

X2 

■ 

A 
I       >  J 

0.O0095H941  -0.0000298211  -0.0000191838 
-0.0000298211  0.00O58U536  -0.0000191838 
-0.0000191838  -0.0000254159     0.0004171769 


11776.1340 

23875.7760 

27264.7400, 


9.96988 ^ 
13.08007 
10.54149 


Deflections   of  the  base   structure:       (anti-symmetrical  con- 
dition) 


mm 
"l*    I     -^ds  =  -31.86760/E 


2       }s^ 


^ds  -  -83.5S46O/2 


mnm-, 
"3"    \     "ll   ds  =  -122.40455/E 


mm       l 
1  1 


11=  \   -Sr^ds  +  J-4  =  59.32237+1018.1785  =  1077.50087/E 


21  AXE 


m  m 
f12*    I     "ir^ds  =  104.96799/E 


/     m^  m, 
f13=    \     -TJpds  *   112.89622/E 


/    mm 

f     -  =      I  21 H  a     -     1 


04.96799/E 
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"glfeda  i   7-1  =  229.23666  +  1610.6911  -  1839-92776/E 


/    ir.-m., 
f23=         TI^8  =  273-73558/E 


f    mm 
f     -    1     -^ids  -  122.39622/E 


|^ds  =  273.73558/E 


■32   i  -ET 

-IJ^is  +  _1  =  419.98676  +  2276.8258  =  2696.91256/E 

>->     I     EjI 


ATE 


put  into  matrix  equation  thus 

'1077.50087  104.96799  112.89622 
104-96799  1839.92776  273-73558 
(  122.89622  273.73558  2696.91256  J 


-h 


-31.86760 
-83.58460 

•122.40455J 


The  solution  from  Computer  output  gives, 
"X 


"1 

X2 
1X3J 


0.0009363737  -0.0000483181  -0.0000342935 
-0.0000483181  0.0005543259  -0.0000342935 
-0.0000342935  -0.0000542412  0.0003777354, 

f  0.0216036 
O.0405956 
[  0.0406099  J 


31 

.86760' 

83 

.58460 

122 

.40455  . 
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("o;   Koneat   diagram  for  anfcicyw.otrical   cable   3tres.;es 
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(a)  Loading  on  the  frar.e  due  to  cable  stresses 
(summation  of  symmetrical  and  anti-symmetrical 
leading  conditions ) 
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Co)  Moment  diagraa  of  the  frame  due  to  cable  stresses. 
Pig.    22 
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(a)  Moment  diagram  of  the  frame  due  to  external  loads 


without  cable  stiffening. 
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tb)  Pinal  Koment  diagr-aii  of  the  frai.ie  due  to  external 
loads  with  cable  stiffening. 

Kg.  23 
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APPLICATION  OP  HAUNCH  MEMBER  THEOREM  TO  CABLE 
STIFFENED  STRUCTURE  ANALYSIS 


In  concerning  the  elastic  behavior  of  a  reinforced  con- 
crete member  at  low  stresses  up  to  about  f'c/2,  the  concrete 
is  seen  to  behave  nearly  elastic,  that  is,  stresses  and 
strains  are  quite  closely  proportional,  because  the  compres- 
sion strain  in  the  concrete  at  any  given  load  is  equal  to  the 
compression  strain  in  the  steel. 

£c  "  fc/£c  =  Es  -  fs/Es   s°  that, 

f«,  "■  •=—  f „  "  nf  ,  where  n  =  E„/E_  is  known  as  the 
0   E   c     c  s'  c 

c 

modular  ratio. 
Let  Ac  =  the  net  area  of  concrete  and  As  =  area  of  reinforcing 
bars.  Assume  a  concrete  beam  is  under  axial  load  P,  that 

P  *  fcAc  +  fsAs  -  fcAc  +  nfcAs  -  fc<Ac  +  nAs>      (33) 

The  term  (Ac  +  nAs)  can  be  interpreted  as  the  area  of  a  fictious 
concrete  cross  section,  called  transformed  area,  which. is 
subject  to  the  particular  concrete  stress  f  result  in  the 
same  axial  load  P  as  the  section  composed  of  both  steel  and 
concrete.   This  transformed  area  consists  of  the  actual  con- 
crete area  plus  n  times  the  area  of  reinforcement. 

If  a  beam  is  under  bending  moment,  it  is  assumed  that 
compression  will  occur  on  the  top  of  the  section,  and  the 
stresses  are  carried  by  concrete  and  compression  steel,  if 
any.   There  is  also  tension  in  the  bottom  of  the  beam  section, 
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and  the  stresses  are  carried  by  steels. 

It  is  assumed  that  the  concrete  is  not  capable  of  re- 
sisting any  tension,  so  the  shape  of  the  transformed  section 
as  shown  in  Fig.  24  (b)  is  acting  by  internal  force  as  the 
stresses  shown  in  Fig.  2k  (c).  Let  At  -  total  transformed 
areas,  and  I.  =  moment  of  inertia  of  the  transformed  area. 
The  maximum  stress  in  concrete  is 


P   Mcr 


"  cmax 


H    xt 
and  the  tension  stress  in  steel  is 

,  d   Me 

At    Jt 
where  P  is  compression  force  of  the  section  and  M  is  the  bend- 
ing moment.   Therefore,  once  the  transformed  section  has  been 
obtained,  the  usual  methods  of  analysis  of  elastic  homogen- 
eous beam  can  be  applied. 
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The  transformed  section  method  is  usually  applied  in  the 
reinfoi'ced  concrete  members.   The  author  wants  to  use  this 
method  in  the  cable  stiffened  structures.   It  is  advantageous 
to  increase  moment  of  inertia  at  the  ends  of  cable  stiffened 
member  without  increasing  concrete  quantities. 

In  Fig.  25  is  shown  a  cable  stiffened  rigid  frame,  the 
cables  on  this  structure  really  work  as  the  reinforcing  steel 
to  the  concrete  beam,  because  the  cables  just  take  part  of 
the  tension  from  the  reinforcing  steel  in  the  beam,  so  that 
the  concrete  and  steel  section  is  reduced  due  to  cable  stif- 
fening.  If  the  cable  sections  are  transformed  into  a  fictious 
concrete  section,  the  beam  is"  really  stiffened  as  a  haunched 
member,  and  this  structure  can  be  analyzed  as  a  haunch  member 
structure. 


Fig.  25 


6k 

It  is  necessary  to  illustrate  some  fundamental  theory  of 
haunch  members  in  order  to  apply  to  the  cable  stiffed  struc- 
tures analysis.   A  haunch  member  structure  can  be  analyzed  by 
general  moment  distribution  method  after  the  stiffness,  car- 
ry-over factor,  and  fixed  end  moment  have  been  determined. 

Now  a  cable  stiffened  structure  as  shown  in  Fig.  25  is 

taken  into  consideration.   This  is  a  long  span  rigid  frame 

with  cables  hanging  on  the  beam  top.   It  is  a  rigid  frame  of 

uniform  section,  span  L,  the  lower  column  height  is  H  and  the 

upper  column  height  is  h.   The  beam  AB  is  hanged  by  six 

cables  equally  distanced  at  points  E,  F,  G,  H,  I,  J,  and 

three  of  the  cables  on  each  side  are  built  in  the  top  of  the 

upper  columns  at  K  and  L  respectively.   The  moment  of  inertia 

of  the  beam  is  I„   and  column  is  I„ . 
b  U 

Assume  the  beam  AB  is  divided  into  equal  span  L  as  shown 
in  Fig.  25,  the  sections  1-1,  2-2,  .  .  .n-n.   Because  it  is. 
considered  that  both  the  cables  and  the  beam  taking  the  loads 
acting  on  the  beam,  so  the  transformed  areas  on  each  section 
of  the  beam  can  be  divided  into  four  groups  on  each  side  of 
the  frame.   That  is,  the  sections  from  A  to  E  of  the  trans- 
formed areas  consist  of  the  actual  beam  sectional  area  and  of 
the  cables;  the  sections  from  E  to  F  consist  of  beam  area  plus 
two  fictitious  concrete  areas  of  two  cables;  the  sections  from 
F  to  G  consist  of  the  beam  areas  and  one  fictitious  concrete 
area  of  cable;  and  the  sections  from  G  to  K  is  the  beam  cross 
sectional  areas,  because  there  is  no  cable  through  these  sec- 
tions. 
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Let  cable  KE  »  LJ  ■  cable  1,  cable  KF  ■  LI  -  cable  2  and 
.  I  le  KG  -  LH  «■  cable  3,  so  the  transformed  areas  on  each 
sectional  group  will  be, 

from  A  to  E,   An  =  A„  +  n(A  ,  +  A  „  +  A  ,  ) 

t-        SX      S*i      Sj 

where  A  ,  ,  A  0,  A  -  are  the  cross  sectional  areas  of  each 
si'   s^'   s3 

cable. 

From  E  to  F,   A2  =  Ac  +  n(As2  +  ^3). 

Frorr.  F  to  G,    A   =  Ac  -t  n'Ugj). 

From  G  to  H,   A,  =  A  . 

In  calculation  of  moment  of  inertia  for  each  section,  it  is 
necessary  to  locate  the  axis  of  center  of  gravity  of  those 
sections  because  the  positions  of  the  cables  on  each  section 
are  varied  and  the  distances  between  the  central  axis  of  the 
cables  and  beam  of  each  section  are  varied,  and  proportional 
to  the  slope  of  the  cables.   For  the  sake  of  simplicity, 
the  central  axis  of  beam  AE  is  taken  as  the  central  axis  of 
each  transformed  section,  because  the  section  of  the  cables 
is  so  small  in  comparison  to  the  concrete  beam  section,  so 
that  the  error  in  location  of  the  central  of  gravity  of  the 
transformed  sections  is  small  and  neglegible. 

The  moment  of  inertia  of  the  transformed  sections  will 
be, 

Tx  "  XB  +  nA8hx  (34) 


fit 

where  the  subscript  x  is  the  location  of  the  section  number 

counted  from  A.   For  example,  in  section  1-1,  the  moment  of 

o       2  2 

inertia  of  this  section  is  I,  ■  ID  +  n(A  =  ih^  ,+A„0hn  +A_,h,  .  )  , 

ID       BJ.  11   S£  12   °-3  13 


where  h. 0  means  the  vertical  distance  of  cable  1  to  the  beam 
axis  at  section  1,  and  so  do  h.  _  and  h   .   For  further  ex- 
amples, the  moment  of  inertia  of  section  5-5  is  Ir  =  I„  + 

n(A  ~h2+A  ,h2  )  and  in  section  10-10  that  Inn  •»  I-  + 
^"-  5 2  s.5  ^^  lu    y 

2 

n(As-jh-,Q-, ) .  Now  the  location  hx  for  each  cable  at  each  sec- 
tion of  the  frame  shown  in  Fig.  25  can  be  written  in  matrix 
form  as , 


N 


111 

M 

'4 

8 

12' 

h/4  ' 

3/4  7/8  11/12 

h 

3 

7 

11 

h/8 

2/4  6/8  10/12 

LI 

2 

6 

10 

h/12. 

1/4  5/8  9/12 

1 

5 

9 

0  4/8  8/12 

0 

4 

8 

0  3/8  7/12 

= 

0 

3 

7 

0  2/8  6/12 

0 

2 

6 

0  1/8  5/12 

0 

1 

5 

0   0   4/12 

0 

0 

4 

0   0   3/12 

0 

0 

3 

C   0   2/12 

0 

0 

2 

0   0   1/12  , 

0 

0 

1  , 

The  moment  of  inertia  for  each  section  can  be  found  by  equa- 
fci<  .-.  (34)  and  the  variable  moment  of  inertias  can  be  expressed 
as 
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V 

f     ,1 

r  nAsl  , 

Xb] 

■  k/ 

nAs2 

+ 

(35) 

-n 

I 

>  nAs3  J 

iJB, 

A  diagram  can  be  ploted  for  tho  moment  of  inertia  in 
each  section  as  shown  in  Fig.  36.   It  is  obvious  that  the 
moments  of  inertia  near  the  supports  are  appreciably  larger 
than  the  central  portion  of  the  beam,  which  has  the  same 
properties  as  the  ordinary  haunch  members. 


2  3  4  5b   i  6  9  iov,  121  jiui  516171 01  §021 22232425202720 

boars   s^ci-ions    (ft.) 


?ig.    26  I      diagram 
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The  bending  moment  in  a  rigid  frame  with  variable  moment 
of  inertia  will  undergo  the  following  change: 

(1)  the  negative  moment  at  the  end  AB  becomes  larger  than 
for  beam  with  constant  I. 

(2)  the  length  of  the  region  subjected  to  negative 
moments  becomes  larger. 

(3)  the  positive  bending  moment  in  the  center  becomes 
smaller  than  for  the  beam  with  constant  I. 

(A  )  the  increase  of  the  moment  of  inertia  at  the  fixed 
end  has  the  effect  of  putting  more  bending  moment 
near  the  supports. 
From  the  characteristics  of  a  haunch  member  as  discribed 
above,  for  the  cable  stiffened  structure,  the  increased  nega- 
tive moment  produced  at  each  end  shall  be  taken  by  the  cables, 
because  the  tension  forces  at  each  cable  attachment  form  a 
moment  couple  with  respect  to  the  ends  A  and  B.   In  other 
words,  the  fixed  end  moments  of  beam  AB  is  really  reduced  due 
to  cable  stiffening,  but  in  this  case  the  bending  moment  of 
columns  at  points  A  and  B  shall  be  increased  due  to  the  in- 
creased tension  forces  of  cables  at  the  top  of  the  columns. 
If  this  structure  can  be  modified  by  extending  cables  from 
points  K  and  L  to  the  ground  E  and  F  as  shown  in  Fig.  27,  then 
both  the  moments  of  beam  and  columns  can  be  reduced  by  cable 
stiffening.   This  kind  of  arrangement  is  better  than  the  one 
shown  in  Fig.  25  unless  it  is  architecturally  displeasing. 
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Pig.  27 


In  order  to  show  the  analysis  of  cable  stiffened  rigid 
frame  as  a  general  rigid  frame  with  variable  moment  of 
inertia,  it  is  necessary  to  derive  some  formulas  for  beam 
with  variable  moments  of  inertia. 

(1)  Fixed  end  moment: 

(a)  For  concentrated  load  condition  general  expressions 
for  the  fixed  end  moments  of  this  kind  beam  as  shown  in 
Fig.  25  can  be  found  by  stating  that  the  changes  in  slope  and 
vertical  deflection  between  the  ends  of  the  beam  must  be  zero, 

that  is 
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(L  Mdx 

Jo  EI 


=  0 


Mxdx 


0= 


0 


Fig.  28 


Where  both  M  and  I  are  function  of  x.   The  bending  moment  of 
this  beam  can  be  expressed  in  two  parts,  that  is,  the  portion 
due  to  the  concentrate  load  acting  on  a  simply  supported 
beam,  and  the  part  caused  by  two  end  moments.   So  the  inte- 
gration across  the  team  must  be  made  in  two  steps  due  to 
the  discontinuity  of  the  moment  curve  for  the  simply  sup- 
ported beam. 
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fedix   H     |LPax(L-x)dx   _    |L^a    _   X^.jyXdX},.   &         (36b) 

Equation   (36a)   can  be  written  in  the  following  form  by  adding 

and  substracting    I     2Li>L]X)   thus 
Jo     LI 

Pb   1  xdx  _   P  (    ,(x-a)d     M     (    dx   _  'A"   D  (   xdx  ,J7    j 

E   ),   I  )Q     LI  MA  )Q   I  L       )0  I  l37a' 

Likewise   the   equation    (36b)    can  be   expressed   thus: 


_  P  W,  =  „  (L2dx  .  %^B  (Vdx  (3?b) 


(^f=^dx  =  ma  r 


L  ,  I      1   I   "A   "A  1   I      L   \   I 

;o  )o  io  o 


Eliminating  *"«-Mg  from  equation  (37a)  and  (37b)  and  solving 

for  M  gives, 


/\dx  (Lx(x-a)H.,.   (Lx2dx  /L(x-s 
l1  )a  X   ^'lo1   )o  1 


'-dx 


h  "  7WTT         /l  ;   2 <3*> 

1  gx  1  xcdx  _  (  I  xdx. 


o1  o1      o1 


A  similar  expression  for  Mg  can  be  obtained  by  changing  a  to  b. 

(b)  For  a  uniform  load  condition;  in  case  a  load  w  is 
uniformly  distributed  over  the  entire  span  of  a  beam  with 
variable  moment  of  inertia,  the  fixed  end  moment  at  the  left 
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end  can  be  found  in  a  similar  manner,  thus 


W, 


I 


xdx  {    x2dx 


0 


rlLx2dxi2 
10 


Lx2dx 

o    x 


Lxdx  - 2 


(39) 


u0*n 


Mg  can  be  found  from  equation  (39)  by  performing  the  inte- 
gration from  the  opposite  direction. 

(2)  Carry-over  factor: 

The  carry-over  factor  is  the  ratio  of  the  moment  induced 
at  the  fixed  end  of  a  member  to  the  moment  causing  rotation 
at  the  other  simply  supported  end.   In  Fig.  29  is  shown  a 
beam  with  variable  moment  of  inertia  subjected  to  the  action 
of  a  moment  M  at  the  simply  supported  end  A  and  resisted  by 
a  moment  of  magnitude  CM  at  the  fixed  end  B.  When  deflection 
at  A  is  zero,  then 


Pig.  29 
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(l-x/L)M  -  £<JM 
0 EI xdx  =  °  <^°) 


Integrating  and  solving  for  C  gives  the  carry-over  factor 
from  A  to  B  is 

r  i  xdx  _  I  x^dx 


Jo.1      h1 

-AB  -    ft^T  (U) 


The  carry-over  factor  from  B  to  A  can  be  found  by  equating 
the  deflection  at  B  to  zero,  as  A  is  fixed,  which  gives 


C3A  -  ~4r ^r, -m rr—  (w) 


(3)   Stiffness: 

The  stiffness  factor  at  one  end  of  a  beam  is  the  moment 
required  to  produce  unit  rotation  at  one  end,  which  is  assumed 
simply  supported  while  the  other  end  is  fixed.   According  to 
the  principle  of  moment  area, 

LfL(i  -  S(i  +  c)]dx  -  e  (1+3) 

o 

when  0  «  1,  K  may  be  substituted  for  M  and  solve  equation 
(43 )  £nd  gives 
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'dx 


1+C  (LxAx 


(44) 


If  the  value  of  the  carry-over  factor  C  in  equation  (41)  is 

substituted  in  equation  (44),  the  expression  becomes 
'Ijx^dx 


K 


AB 


(45) 


Similarly,  the  stiffness  Kda  can  be  found  by  substituting 
Cg^  to  equation  (44)  and  gives 


VBA 


-JM 


jo  I  )o~ 


(46) 


(4)  Lateral  stiffness: 


Pig.  30 


The  lateral  stiffness 
is  the  moment  produced  at 
one  end  of  a  member  fixed 
against  rotation,  if  the 
other,  also  fixed  against 
rotation,  is  displaced 
laterally  with  respect  to 
the  first  end.   In  Fig.  30 
is  shown  a  beam  with  variable 
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moment  of  inertia,  in  which  the  stiffness  K,_,  K,,,  and  carry- 

AB   BA 

over  factor  C,_,,  C„.  have  already  been  determined.   The  end 

A  B '   BA  J 

A  of  the  beam  has  been  displaced  an  amount  I\    .   It  is 
desired  to  find  fixed  end  moments  K  and  M  .   Because  there 
are  no  forces  acting  on  the  beam,  so  the  moment  at  a  point, 
at  a  distance  x  from  the  end,  is 

Mx  "  MA  "  I{MA  +  V  M) 


According' to  the  principle  of  moment  area,  the  parallel 
tangents  to  the  elastic  curve  at  the  end  can  be  expressed 

as 

)o(Ma  '  i(ma  +MBn  nr  =  °  (48) 

The  vertical  displacement   -A  between  the  ends  can  be  expres- 
sea  as 


i 


L 
)Q(\  -   |(MA  +  MB)  xdx  >  _A  U9, 

Equation  (48)  and  equation  (4  9)  can  be  written  thus: 

KAL  pf  "  (MA  +  MB)   jLx^=  0  (50) 

MAL  (L^  -  (MA  +  MB)  (Vdx  .  _EA  (51) 


IK 


Solving  for  M,    gives, 


xdx 


M,    -EA-7J tyr-j 71 5-  (52) 

A  /^dv    f  S^d*        ,  <  Lxdx   -,  * 

lo  *    J 


pp^ 


and  also   can  he  written: 


K 


Lx2dx  (Lxdx 


K,   -A-* fe4 Tr 5  x     i?  ,.  (53) 


Fto   fLx2dx t  pxdxl2  (Lx2dx 

)o"r)o"T"  l)oiJ    jo1 


where  the  first  fraction  in  equation  (53)  is  the  stiffness 
and  the  last  fraction  is  equal  to  (1+C)*  1,  therefore 

\  '  I   W1  +   CAB}  (54) 

Similarly 

MB=IKEA(1+CBA)  (55) 

The  above  derived  formulas  are  the  necessary  fcators  for 
solving  frames  with  variable  moment  of  inertia.   A  numerical 
example  will  be  given  to  describe  the  method  of  applying 
haunch  member  theories  to  solve  cable  stiffened  structures. 
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NUMERICAL  EXAMPLE  3 

Find  the  fixed  end  moment  of  the  cable  stiffened  rigid 
frame  as  shown  in  example  1.   Consider  the  horizontal  cable 
stiffened  girder  as  a  haunched  member.   The  frame  is  uniform 
loaded  with  2.25  k/ft. 


5  1 6  I  7  1 8|  910lll2l314l5'l6l7l819202122232425262728 


Ib=20.833 
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Fig.  31 


In  the  frame  as  shown  in  Fig.  31.  the  horizontal  cable 
stiffened  girder  is  divided  into  28  equal  intervals,  as  shown 
in  the  figure,  with  sections  from  1-1  to  28-28,  each  inter- 
val being  5  feet  long,  the  frame  is  2'wide. 

The  first  step  is  to  locate  the  new  axis  of  each  section 
for  the  transformed  areas.   For  the  uniform  horizontal  girder, 
the  central  axes  are  all  the  same  throughout  the  sections, 
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and  are  located  at  the  half  way  of  girder  depth.   Duo  to  the 
additional  transform  area  of  cables,  the  central  axis  for 
each  section  will  move  a  little  amount  upward  from  the  central 
axis  of  the  girder.   The  eccentricity  for  locating  the  new 
axis  will  show  in  the  following  table. 

The  second  step  is  to  calculate  the  moment  of  inertia 
in  respect  to  the  new  axis  for  each  section.  It  is  simple 
Without  mathmatic  expression. 

The  third  step  is  to  calculate  the  carry-over  factor  and 
stiffness.   First,  a  unit  moment=l  is  applied  at  section  1, 
then  at  each  section,  and  the  corresponding  moment  due  to 
the  unit  moment  can  be  found  by  proportion,  so  that  the 
diagram  for  moment  divided  by  moment  of  inertia,  (M/l)  can 
be  plotted  as  shown  in  Fig.  32. 

The  slopes  of  the  M/l  loading  diagram  can  be  found  by 
summation  of  those  ordinates  from  section  1  to  23,  and  take 
moment  in  respect  to  the  left  end  of  section  1.   Let  the 
slope  9q  be  the  reaction  of  the  girder  loaded  with  K/l  dia- 
gram at  point  0,  and  62g  be  the  reaction  of  the  girder,  for 
the  same  loading  condition,  at  end  point  28,  and  get  the 
relationship  between  stiffness  S  and  carry-over  factor  C  as: 

SxG0  -  C£x92g  =  1 
Sx62g  -  CSx90  =  0 

The  results  of  the  above  simultaneous  equations  gives  the 
carry-over  factor  and  stiffness.   From  the  computer  answer, 


•  ' 
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Section 

Transform 

rABLE  1 

Essentricity 

Moment  of 

No. 

Section  Areas 

Inertias 

(ft2) 

(ft) 

(ft*) 

28 

10.75006 

1.39545270 

299.923540 

27 

10.75006 

1.13225350 

222.664900 

26 

10.75006 

0.96906436 

161.444290 

25 

10.75006 

0.75537025 

116.261710 

24 

10.75006 

0.54267611 

87.117154 

23 

10.50004  - 

0.45633405 

66.740130 

22 

10.50004 

0.35717016 

50.746025 

21 

10.50004 

0.25795627 

39.059457 

20 

10.50004 

0.15874237 

31.680424 

19 

10.25002 

0.12196035 

26.931048 

13 

10.25002 

0.08130727 

23.543246 

17 

10.25002 

0.04065369 

21.510563 

16 

10.00000 

0.00 

20.833333 

15 

10.00000 

0.00 

20.833333 

14 

10.00000 

0.00 

20.833333 

13 

10.00000 

0.00 

20.833333 

12 

10.25002 

0.04065369 

21.510563 

11 

10.25002 

0.08130727 

23.543246 

10 

10.25002 

0.12196085 

26.931048 

9 

10.50004 

0.15874237 

31.680424 

3 

10.50004 

0.25795627 

39.059457 

7 

10.50004 

0.35717016 

50.746025 

6 

10.50004 

0.45638405 

66.740130 

5 

10.75006 

0.54267611 

87.117154 

4 

10.75006 

0.75587025 

116.261710 

3 

10.75006 

0.96906436 

161.444290 

2 

10.75006 

1.18225850 

222.664900 

1 

10.75006 

1.39545270 

299.923540 
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o.ooo44Z4A 
o.ooo\&c4o 
0.00000000* 
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there  are: 

carry-over  factor  =  0.73432514  stiffness  -  2.03819740 

Because  the  frame  is  symmetrical,  so  the  carry-over  factors 

are  same  for  both  ends . 

end 

The  fourth  step  is  to  calculate  the  f ixedAmoment  in- 
fluence line.   This  is  done  by  the  conjugate  beam  method. 
First  cut  the  fixed  end  0  and  put  a  moment  and  a  reaction 
there.   This  cantilever  beam  is  loaded  with  a  moment  and  an 
upward  force  at  point  0,  so  the  M/l   loading  diagram  for  this 
cantrilever  beam  can  be  found  as  described  in  step  2.  The 
reaction  of  the  fixed  end  for  the  M/l  loading  is  the  summa- 
tion of  the  loading  diagram.  So  the  slope  theta  is  the 
loading  diagram,  as  shown  in  Fig.  33- 

The  deflection  of  the  cantilever  is  the  fixed  end  moments 
for  the  M/l  loadings.   The  parts  of  this  deflection  curve 
enclosed  by  a  straight  line  through  both  ends  of  the  curve 
is  the  fixed  end  moment  influence  line  of  the  girder,  which 
is  the  same  for  both  ends  due  to  symmetrical  condition.   The 
deflection  of  the  cantilever  and  the  fixed  end  moment  influ- 
ence line  of  the  girder  are  shown  in  Fig.  34. 

The  fifth  step  is  to  calculate  the  fixed  end  moment  due 
to  external  loads,  which  is  the  integration  of  the  moment 
influence  areas  times  the  uniform  loads,  and  it  gives  FEM  = 
3996.9413260875  kips-ft. 
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S>-Z013Z78G 

^■4ZZs>&&32> 

^■54Z54W3 
\p.£,ZZ833Z3 
^■TZZ^SoZt 

^p.?T33oZ&Z- 

\l.  i&033Y3O 
iZZ&M^O 

j.44y4I$oO 
\J<49267?70 

\l.53l04400 
\1.5G383010 
15359I230 
\i.6Zo3&<£>50 
\/.£>Z<2>67010 

\l.&&l  l&l4o 


' 
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0.OOOO449O 
\4.l23l?l40 

3.05397950 
U.7J742400 

IP.  OS>6o<350 
J3.oT  170400 

P0.&4&75&0O 
\Z2.e5o2o3oo 
\Z4. 12055Z0O 
\24.6&54IOO 
&5024I000 
Z453447?>oo 
Z3.47G76T°° 
ZI.&46 19300 
Z0.053&ZZ00 
1T9I03Q&00 
1 5.6071 0500 
1327153500 
1Q97I205OO 


~  3.394 13300 

\^-9a9790lO 

J5-4I5&4740 
\4.05zQi1 10 

1 2.99044530 
!2.looo&430 
'I-4I3&IOZO 
'0.34244437 
fe-39135,25^ 
o,  0000  <S  1 1 r_5 
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The  final  step  is  to  balance  the  fixed  end  moment  by 
moment  distribution,  which  gives: 


0.0)1,9  o .  -K  i 


+3391)..  0032 

+339i|.0032 


-3996.9:i.133 
-I-6C2.  93&1U 
-339^.0032 


K'=2. 0382x0.5=1 .0191 


«f 


+1697.0016 

(HOTS:  These  moments  are  in  ft-kips) 

The  fixed  end  moments  are  somewhat  larger  than  the  foregoing 
examples,  because  the  negative  moment  for  a  frame  with  vari- 
able moment  of  inertias  is  larger  than  for  frame  with  constant 
moment  of  inertia.   These  results  are  the  approximate  answer, 
because  this  is  an  approximate  analysis.   The  more  smaller 
interval  sections  are  divided,  the  more  accurate  the  answer 
it  gives.   The  moment  of  the  girder  end  is  the  total  moment^ 
in  designing  the  horizontal  girder,  the  reinforcement  for 
this  moment,  deduct  the  cable  areas  is  the  actual  steel  areas 
for  the  horizontal  girder  reinforcement  because  the  cables 
take  part  of  reinforcing  function  from  the  girder  reinforce- 
ment . 
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This  method  is  not  suggested  for  actual  design,  because 
It  is  qui  to  laborious  to  calculate  section  by  section,  and 
it  is  impossible  to  give  high  accuracy.   The  purpose  of  this 
example  is  to  introduce  one  kind  of  method  for  cable  stif- 
fened structure  analysis. 
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CONCLUSION 


A  cable  stiffened  structure  usually  is  a  composited 
structure.   The  cable  can  be  used  as  stiff ener  for  any  kind 
of  reinforced  concrete  rigid  frame,  steel  frame  or  truss  as 
well  as  timber  structures.   This  kind  of  composite  structure 
always  combines  different  stress  conditions  in  one  structure: 
the  cable  takes  tension  only;  the  frame  can  take  moment,  shear 
and  thrust;  and  the  truss  is  assumed  to  take  direct  force 
and  shear.   Therefore,  a  cable  stiffened  structure  always 
has  a  high  degree  of  indeterminacy.   In  this  condition,  the 
best  way  is  to  use  the  principle  of  superposition  to  solve 
the  simultaneous  equation. 

In  this  report  the  author  presents  a  typical  long 
spanned  rigid  frame  with  six  cables  hanging  on  the  top,  which 
has  been  analyzed  by  both  superposition  method  and  haunched 
member  method.   Checking  the  numerical  answers  from  the 
numerical  examples  shows  that  the  fixed  end  moments  are  some- 
what different  due  to  different  approaches,  but  the  method 
of  superposition  is  more  accurate  than  the  haunched  member 
method  because  the  haunched  member  is  calculated  by  taking 
small  sections.   If  more  smaller  sections  are  divided,  more 
accurate  results  can  be  obtained. 

Three  numerical  examples  have  been  presented  in  this  re- 
port.  The  first  one  is  calculated  by  superposition  method 
and  the  frame  is  considered  to  be  loaded  with  vertical  dead 
and  live  loads  only.   The  frame  is  analyzed  in  symmetrical 
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condition  and  has  a  moment  reduced  by  the  cable  about  one- 
fourth  of  the  original.   Here  the  author  has  tried  to  con- 
trol the  fixed  end  moment  to  bo  reduced  in  an  amount  as 
desired.   It  is  done  by  prestressing  the  cables  in  a  definite 
vertical  lengthening.   The  numerical  results  for  this  pre- 
tension work  are  satisfied.   The  second  example  shows  the 
frame  under  lateral  wind  stresses,  assumes  no  live  load  on 
the  top,  and  analyzes  it  by  superposing  the  symmetrical  and 
antisymmetrical  conditions.   The  results  show  that  under 
normal  dead  load  condition,  the. cables  are  still  all  in  ten- 
sion due  to  lateral  wind  stresses  and  the  frame  is  little 
influenced  by  the  live  load,  because  the  span  is  very  large 
and  the  dead  load  of  the  structure  is  much  larger  than  the 
live  load.   The  third  example  shows  that  the  fixed  end  moment 
of  the  frame  is  calculated  by  assuming  the  horizontal  girder 
is  haunched  by  cables,  which  gives  the  fixed  end  moment  for 
haunch.   After  moment  distribution,  it  is  found  that  the 
girder  end  moment  is  larger  and  the  column  moment  is  smaller 
than  in  the  first  example. 

This  report  has  presented  computer  method  for  solving 
problems.   The  major  part  of  the  three  examples  are  calcu- 
lated by  computer.   Those  examples  also  can  be  calculated  by 
hand,  but  with  much  labor.   It  is  evident  that  it  is  more 
convenient  to  solve  simultaneous  equations  by  computer  in 
analyzing  high  degree  inderterminate  cable  stiffened  structures. 
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In  conclusion,  the  cable  stiffened  structure  is  very 
useful  in  modern  architectural  structure  design.  There  are 
Biany  possible  way's  of  structural  analysis  for  cable  stiffened 
structures,   it  is  the  author's  opinion  that  each  method  in 
a  different  approach  can  get  agreement  in  result,  but  the 
time  consumed  and  percentage  of  accuracy  should  be  taken 
into  consideration.   Most  of  the  calculations  in  this  study 
have  been  by  computer  rather  than  slide  rule. 
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APPENDIX  A 
NOTATIONS  AND  ABBREVIATIONS 
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...'  MMX  A   NO'i'ATI'    AND  i   I'-isSVIAT-tOMS 

A 

c 

Cross  soctioal  area  of  concrete  girder 

A  . 

Cross  sectional  area  of  steel  or  cable 

A3 

External  rigidity  of  cable 

Jb 

Moment  of  inertia  of  beam 

I 

0 

Moment  of  inertia  of  column 

T^ 

Moment  of  inertia  of  the  transformed  area 

SI 

The  flexural  rigidity  of  concrete  bean 

>t 

coefficient  for  extensional  rigidity  of  the 

effective  straight  cable 

bT 

The  increment  of  cable  tension 

£a 

A  unit  deflection  at  point  A 

A! 

Elongation  of  cable 

K 

Stiffness  of  the  frame 

K' 

-he  modified  stiffness 

L 

The  actual  length  of  cable 

X 

The  cable  stress 

u 

The  strain  energy 

■( 

Moment  distribution  due  to  the  action  of  external 

load  on  the  base  structure 

m  X 

n  n 

The  moment  distribution  due  to  the  action  of  the 

arbitrary  constant  X  alone 

fij 

The  influence  coefficient  of  deflection  on  base 

structure 

u          : 

The  deflection  on  the  base  structure  due  to  appli- 

ed load 
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[FJ 

:  Flexibility  matrix 

txj 

:  Column  matrix  for  cable  stress 

Luj 

:  Column  matrix  for  cU.3Dlacsment 

[L] 

:  i.owrar  traingular  matrix 

K 

:  Koi-ient  In  kip-ft. 

Ah 

■  'ihe  vertical  elongation  for  cable  prestressing 

CA. 

Carry  over  factor  from  A  to  B 

n  "  VSc 

'.he  Modular  ratio  oi'concrete  and  steel 

nAs 

Transform  area  for  cable 

'  or  ft . 

foot (feet ) 

[!  or  in. 

inch 

k 

kip 

so. 

s  quare 

,r  or  lb 

pound 

0 

Skiameter  of  cable 

e             : 

angio 

St/'         : 

Icip  per  foot 

-, 
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COMPUTER    PROGRAM    .    DATA    AND    OUTPUT 
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PPC 

r 


1 
C 

I 

i  = 
ci  = 

DC 

RFA 

REA 

D 

REA 

Al  = 

A2  = 

RcA 

DC 

Z(  I 

DC 

J=I 

ZW< 

A=Z 

A=A 

A  =  A 

A=A 

PUN 

CCN 

CON 

FOR 

STG 

PN'D 


1 


] 
2  I 


FOR  FI 
M  DEFLE 
EMS  I  CM 
. 
. 
201  LL= 
D,XA,X3 
D.UI  I) 
100  K=l 
D.YA.Yd 
25.*tXA 
20.»XBC 
D.I Y( I ) 
10  1=1. 
)=X( I )* 
2(1  1  =  1, 
+  1 
I  )  =  IX(  I 

ii)+2.» 

'--]  ./J. 

+5.*(Z( 
/BI+IA1 
CH  30, K 
T1NUE 
TINI  I 

"ATI  I5.| 
P 


RST  EXAMPLE 

CTICN  AT  CABLE  PTS.  DUE  TC 
XI15)  ,Y(  lb)  ,Z  I  15  )  ,/."(  15  I 


VERTICAL  I 


JIN-SIEN  GUC 


A 

1  ,5) 

bC 

A+XB*YB+(XA+Xb)*(YA+Yb) )/3. 

BC/3. 

=  1.-5) 


I  ) 


)+X( J) )*( Y( I )+Y( J) )/4. 
(Z(2)+Z.<  3)  )+4.*(Z!  ■( 1  )+Z  (2 )+Z- 

41+ZI5) ) 
+A2J/CI 

,A 


F18.8) 


I  3)  1+ZI4) 
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DATA 

11813.61250       -3627.22500       0.001 

-3627.2250      -927.2250      +672.77' 


+    .99702 
+12.14797 
+2.52I  87 

+12.83726 
+  /i  .'  1246 
+  i  0.94649 
+    .997    ? 
+12.14797 
+    .99702 
+12.14797 
+2.5208? 
+12.83726 
+4.     1246 

*  ]    1.94649 

♦  2.521  87 
+1 2.83726 
+?.52    87 
+1 2.8?776 
+4.     1 246 

+  i  3.94649 
+    .99702 
+12.14797 
+  4.1    1246 
<  ]     .   14649 
+  '. .     1246 
I  I     .  '=4649 

.  2087 
+12.83726 
+  .99702 
+12.14707 
OUTPUT 
C   C   BEAM 


-  .     1403 
-1 .9  14 

-5.    4174 

+  3 .  8  9  1  3  9 
-8.024<  1 

+4.6227: 
-1.99 

-1.994  ! 
-1.994 

-1.99403 
-5.04174 

-8.02491 

+4.6227 

-  .     -174 
+3.89139 

-5.    4  174 

+3.8'    ' 
-8.02491 

+4.62270 
-1.99403 

-1.9 
-8.J2491 

+4.62270 
•     .'491 

+4.6227l 

-  .      't  1  74 
+3. 

-1.994 
- 1 .  9 94<  3 


+1772.7750   +1885.2750 


-14.14160 

-) .09403 
-17.8 

-!   .     4174 
-18.971  41 

-1.7 
-14.14160 

-1.9 
-14.14160 

-1.99403 
-17.87 

-:   .-4174 
-10.  9  M4i 

-1.70128 
-17.87 

-5.    4174 
-17.8 ( 

-'  .    4174 
-18.9714 

-1.70128 
-14.14160 

-1.99403 
-18.97140 

-1 .701^b 
-18.97140 

-1.70128 
-17.87 

-5.1  4174 
-14.14161 

-1 .99403 


- ' .      «403 


-     .    4174 


-8.024    i 


-1.9'   ' 


-5.04174 


-     .     2491 


-1.99403      -1.9 


•1.99403 


-5.041  /4 


-8.C  249] 


-5  .'4174 


-5.04] 74 


-8.     2491 


■1  .'.-9403 


-8.02491 


-8.r.2491 


-5.    4  174 


J  9  4  0  3 


-5.    -.174 


-8. 


-     .     4  174 


-l  .4174 


-8.02491 


■1.99403 


-b.l,<:491 


-8.02491 


,1  4174 


-1.99403       -1.O9403 


LEGTICN    AT    C«qLE    PTS.    DUE    TC    VERTICAL    LOADS         JIN-SIU 


1  -1807    .64  j       . 

2  -31524. 613oOI  I 

3  -31746.20100000 

1  55. ..8971900 

2  72.7558470' 
66.19554  1   II 

125.375': 
123.314/' 
72.7^5F,i7r. 

1  134.5656! 

2  123.314/'  - 

3  66.19554 

STCP   END  OF  PROGRAM 


AT   STATE  I   i    30  +01  LINES 
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PROGRAM  2  FOR  SECOND  FXA^PLE 

C   r    I"  DEFLECTION  AT  CABLE  PTS  DUE  TO  WIND  LOADS 

DIMENSION  XI  15.4)  »Y(  16.4)  ,0(4,4)  ,  PX(  15)  ,PYC16) 
1=20.833 

CI=288.*- 

DO  ZOO  LL=1»2 

IFILL-1 111,11,12 

11  PUNCH  61 
GO  TO  3  1 

12  PUNCH  82 

31  DO  100  N=l,4 

RFAD,  W.TY 

?FADf  Y< l,N) ,Y( Ll.N) 

READ,  X ( 1 , N ) »X ( 3 »N ) , X ( 5 »N ) »X ( 7 , N ) , X ( 8 , N ) 

l-'  =  .( 

Df=(Y(  11  ,N|-Y(  1  ,N)  1/K  . 

DC  10  1=2,11 

H=H+5. 

J  =  I  -1 
1   Y(  I  ,N)=Y(  J,N)+DY+(50.-H)*<v*H/2. 

H«0.0 

ST=TY/2o. 

DO  2r  [  =  1,5 

J=17-I 

Y  (  J,N  )=-«'*H*«2/2.  +  ?T*H 
20  H=h+5. 

DO  3C  1=2,6,2 

J  =  I-1 

<=I+1 
3   X  (  1  ,M  =  CX(  J,N)+XU»N)  )/2. 

(  Oft  1  [NUE 
:  30C  IT=1,4 

IF(IT-lT310.310i 
510  r ' !  =  I  T  + 1 

GO  TO  330 
32C  IN=IT 
33  [D=IN,4 

DO  40  1=1,16 
4r  PY( I  )=Y(  I  ,IT)*Y(  I , ID) 

DO  45  1=1,8 
45  PX(  I  )=X(  I , IT)*X(  I ,  ID) 

Al=(  . 

DO  50  1=2,1  ,: 

J=I-1 

K  =  l+] 

5C  A1=A1+PY( J)+4.*PY( I 1+PYIKI 


JIM-SI  EN  GUO 
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DO  6U  1=13.15.2 

J  =  I  - 1 

<=!+l 
60  A1=A1+PY( J )+4.*PY< I )+PY (K.) 

Al=Al*5./( 3.*CI 1 

A?=  . 

DC  70  1=2,6.2 

J  =  I  -1 

K=I+1 
7'   ;.  ?  =  A?  +  PXI  I  ]+4.*PX(  !  1+PXIK) 

A2«(A2*10./3.+(PX(7)+PX<8) >*5.)/BI 
D( IT»ID>=A1+A2 
3«       PUNCH   8    »IT>ID«Ot  ITflD) 

3£       CONTINUE 

CONTINUE 
8C  FCRMATI2I5.11HDI  FLEC1  I0N=,K16.o) 
PI  FORMAT  (19HSYMMETRI.CAL  LOADING) 
P2  FCRMAT(23HANTIStfMMETRICAL  LOADING) 

STOP 

DATA 

+  .1       -?  . 

+161 i .68875   -3244,2] 

-3224.21083    -844..Z1083    +755.76917    1555.78917    1675.78917 

-14.14160 
+  .99702   -1.99403 
+12.14797   -1*99403   -1.99403   -1.99403   -1.99403 

-17.879 
+2. 52087   -5.04174 
'  .  :-726   +3.89139   -5.04174   -5.04174   -5.04174 

-18.97141, 
+4.1   i  24S       8.    2491 

M     .     4649      +4.6227-       -1.70128      -8.0249J       -8.02491 
+    . 1  -20. 

-91.23976       -20.03165 
-    •■  ^16c      -0.02366      -0.01577      -0.    0789       0. 

-14.14161 
+2.63763  -5.27525 
+=.86675      -7.8081H       -4.68486      -1.56162       0. 

-17.879^0 
+  ■^.96718  -11.93435 
+5.94465      -4.6929]       -15.3313(       -5.10992 

-18.97141 
+o.^B(^q  -16.76065 
+2.21075      -4.184P6       -11.7001]       -18.65557  .00 


'   ,T  =  '1 

c  r 

DEFLECTION 

AT  CAbLF  PTC-  DU 

TRICAL  LOAPINC- 

1 

2DEFLECTICN* 

-1177  6.  1  -  ■• 

1 

3DEFLECTI'CN» 

-23875.77c 

1 

4DEFLECTICN* 

-27264.74  ;  Ut   0.1 

2 

2DEFLECTICN* 

35.60181 

2 

3DEFLECTICN= 

56.02361900 

2 

4DEFLECTICN* 

3  !  .   /' 

3 

3DEFLECTICN* 

1  7.8320510c 

3 

^DEFLECTION" 

i   /.2747PC0C 

4 

^DEFLECTION" 

129.05937000 

ftfjTISYMMETRICAt  I  CADII 

i 

2nFFLrrTir-'- 

-31.86760300 

1 

3DEFLECTICN= 

->83.  58460200 

1 

:  LrCTION= 

-122.40455000 

2 

2HFFLECTICN= 

59.32236700 

2 

3DEFLECTI0Na 

lc4. 96799000 

2 

4DEFLECTI0N= 

112.8962200c 

3 

3DEFLECTICN= 

229.2366' .:  0( 

3 

4DEFLECTICN" 

273.73558000 

4 

4  DEFLECTION 

419.98676000 

96 


WIND  LOADS       JIN-SIr' 


STOP   END  OF  PROGRAV  at   STATEMENT  0(  82  +  01  LINES 
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3  FCR  FIRST  EXA'-VL! 
C   C   CALCULATION  CF  CABLE  STRESSES  BY  IATRIX  INVERSION     JIN-SIEN  GUC 
DIMENSION  A(20r20)*B<39),C(39) 

:  n,N 
MN=N-3 
]   READ. I IA( I»J),J*}»N> ,I=1»N) 
A<l»l>=l./Atl,l) 
DC  110  M«l »MN 

<='-:  +  1 

5C  DC  6l  I  =  1 »  -i 

D( 1 )=0 

DC  6>.  J=]  iM 
60  B( I )=B( 1 ) +A( ] >J)*A ( J»<) 

D=0 

DC  70  I  =  1,m 
70  D=D+A(K» I)*BI I! 

D=-D+A(K,< ) 

A(K.K)=1./D 

do  er  i=i, m 

A( I ,K)=-b( 1)*A(K,K) 

DC  90  J=1,M 

C(J)=0 

DC  9C  1  =  1  ,'■ 

C( J)=C( J) +A1K.I )*Al I ,J) 

DC  100  J=1,M 

I  A (K ,J)=-C( J)*A(K,K ) 
DC  110  I  =  1,M 

PC  110  J=1,M 

II  1  (  I  »j)=A(  T  ,J)-R(  I  )*A  (K»J) 

PUNCH  67 

DO  170  1=1, N 

DC  170  J=1,N 
17   PUNCH  66,  A( I ,J) 

READ,  <B<  I  I  ,1  =  1, N) 

DC  2.u  1=1, N 

C( 1 1=0.0 

DC  200  J=1,N 
?         C( I )=C( i )+A( I »J)»B (Jl 

PUNCH  68 

PC  300  1  =  1  ,M 
30   PUNCH  66,  CM  I 

66  FCRMATIF?     .10) 

67  FORMAT!  M4X.16MMATRIX    INVERSION] 

68  FCRMAT!/14X,14HCArfLE    STRESSES) 
STOP 

END 
3 

1^7^.268219            72.755847  66.193i47 

72.755BA7       1 736. -'6669  J  123.314750 

66.195547         123.31475C  2411.491481 

31524.61  31746. 20Ko 
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OUTPUT 

r      r      CALCUI   ' T ' "'     OF    CAPLF    STRESSES    BY       ATRIX    [NVERSICN  JIN-SII 

IftfRIX    INVERSION 
.00    935743* 
375273 

?~767? 

375273 

5796196 

-.  237672 

02  '76  72 

'   286i-95 

I  £ 7T65 

CAI  LE    STRESSES 
+  1  a . 9 719331 
+  i  (  .33  >623 

+  11.9       3  1  r- 

STOP   FND  OF  PROGRAM  AT   STATEMENT  0068  +  01  LINES 
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EXAMPLE 

CALCULATION  OF  CABLE  STRESSI     I    rRIX  INVERSION     JIN-SIF 
!    '  :  ""  At  2  ,;■     i  •  (39). CI 

DO  350  LL  =  1  ,2 

IF    ILL-]  I     11  .11  .12 
i l    PUNCH    81 

60  10  1- 
12  PUNCH  62 
1       READiN 

NN=N-1 

""P.I  (A(I  »J)  ,J  =  1  »N)  ,1=1. N) 

A<]»1 1=1. /All. 1) 

DC    11C     «=1  . 

<='•  +  ] 
5i     DO    60    1  =  1  ,  ■■ 

B(  I  )  = 

DO    6'-    J=  1  ,  ■  ■ 
( r >=Bl I )+A( 1 .J)*AI J.K) 

D«C 

DO    70     1=1 »M 
70    rw  +  .'  |     .  '  i  -     (  I  ) 

D=-D+A (K.Kl 

A(K.K)=1./D 

DO  8       I  =  1»M 
pn    A  I  I ,K)=-BI 1 1»A(K.<) 

DO    90    J=l «M 

C(J)=( 

(  b  3  i  M 
9u    C(J)=C(J)  t-A(K,I  )*/>(  I.J) 

DO    1'        J=1»M 

I  A  K.  J)=-C(  J)*»"-  (K.K) 
^0    ] 1        1  =  1  .N 

DO  ,11<     J=l  »M 

II  A  (  I  ,J)=A(  I  ,J)-M  I  )*A  (K,J) 
PUNCH    67 

DO    170     1  =  1. M 
DC    17(     J  =  l  ,N 
i 7       PUNCH    6 6 1     A (  I , J ) 
■  '   .  '     I  I  I  .1  =  1  .I'll 
DC    20C     1  =  1, N 
C  (  I  )=     . 
DO    20C    J=1»N 
?    i     C(I )=C<1)+A( I.JI*Bf J! 
PUNCH    6b 
DO    300     I*1»N 
3i  I     PUNCH    66,    cm 
66    FCRMAT(F2<  ,10) 

ft?  FCRMAT(/14X»16HMATRIX  INVERSION) 
«  FO  MAT  !  /l  4X,  1  4Hf«n|.f  STRESSES) 

»1  FCRMATI33 -  2  FOR  SYMMETRICAL  LOADING) 

P2  FORMAT! 38HPR06LEM  2  FOR  4NTI-S1    I  ICAL  LOADING) 
35   CONTINUE 
STCP 
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)  !  T  » 

3 

]  !   .  '■■            .  2362  51.81 

56.02362    1718.5231!  1(  /.2747K 

51.87        1  '  7.27478  2405.  >B>.   17 

-11776.1340   -2381  i  -27264.: 

3 

+  1  77. -V  87    +104. 96799  +112.89622 

t   't.96799   +183  .  2776  +273.73558 

+  112.89622    +.273.73558  +Z6   . 

-31.8676'.,     -83.58460  -122.40455 
0 1 '  T  P  U  T 

C   C   CALCULATION  OF  CAcLfc  STRESSES  L'l  MATRIX  IN\  JIN-SIEN  ObC 

PPOBLE'-'  2  FOR  SYMMETRICAL  LOADING 

MATRIX  INVERSION 

.9514941 

-.01  Oi  298211 

-.;  i  I  191838 

-.1   •  i  298211 

3844536 

0191838 

-.00   1  -J  183b 

2  34159 

4171769 

i   5SES 
-9.969o61 i 
-13.08C  68 
-1  .541486 
PpCLFI'  2  FOR  ANTI-SYMMETRICAL  LOADING 

MATRIX  INVERSION 
.  9363737 
-.0000483181 
-.1  I  0342935 
-,00u0483181 
3543259 
-.OOC  1^2935 
-.''    342 

r-  4  24  12 
,i  03777354 

'  SES 
-.02160  3646 
-.  4 

-.04U6U9948U 
STOP   I   3  01  PROGRAM  Al   STATEMENT  00.82  +  01  LINES 
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PRC 

c 


10 


20 


GRAM 
C  CA 
C  I 
Dt 
RE 
AS 
DP 
DC 
DC 
P( 

:  c 
P( 

DC 
J  = 
P( 

r 

Do 
J  = 

P( 
p 
DC 
J  = 
P( 
DC 
DC 
PI 

P] 

DC 
P  I 

Bl 

B( 
DC 
Al 

DC 
Al 
DC 
C  ( 
DC 

r  ( 

DC 
DT 
DC 

DC 
P2 
DC 
51 
DC 
SI 


6C 


61 


62 


65 


66 


67 


68 


7 


71 


4  FOR 

LCULA 

I   ' 

AD.H, 
=A$#9 

=  H/4. 

10  I 
It,  J 

I  ,J>  = 
?       J 

1,J)  = 
i,      1 

1-1 

I»1J  = 

=  H/8. 
4C  I 

1-1 

I  >?)  = 

=h/12 
50  I 
1-1 
I  ,3)  = 

60  I 
6u  J 

(I.J) 

61  I 
(  I  .1  ) 

6?  I 
11*21 

65  I 
I  )=AS 
4)=AC 

66  I 
I  )  = ..  . 

66  J 
I  )=A( 

67  I 
I)=0. 

67  J 
I  )=C( 

68  I 
(  I  )=C 

70  I 
70  J 

I  I  .J) 
/]  1 

tK)=C 

11  J 
(K)=S 


THIRD  EXA  "I 
TICN  Ch  HUNCH         DUI  TO  CABLE  STIFFENING     JIN-SIEN  GUC 
M3  )  fBl30)»C(15)»DT(15)  »BM(2  ).  [  (30)  >T(30)  »PI  15»3  ) 

CN  P2( 15,3) .PI ( 15,4) 
AS,DX,Af,CI 


=  2.12 

=  1.3 

• 
=  1,3 
n 
=  2,4 


P( J.ll-DP 

=  2,P 

P( J.21-DP 

=  2,12 

P( J.31-DP 

=  1,12 

=  1,4 

=  1.0 

=  6,12 

=  . 

=  H  .12 

=  ■  • 

=  1  ,3 

1,12 


1,4 

)+Pl( I >J)*B< J) 

1,12 

1  »3 

)+P( I »J)*B(J) 

1,12 

I ) /A  !  I  I 

=  1,12 

=  1,3 

=P ( I , J)**2 

=  1,12 

I 

=  1,3 

I (K I+P2 (K»J)*B<J) 
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DC    72    <=1 ,  12 

7?    SI(K)=?I(K)-MK>*<"IT<<1**?) 


PUNCH  1 

PUNCH  1  1 , (SI (K)  > 

K=1.12) 

PUNCH  H  2 

PUNCH  K  1  ,  (  JT(  I  )  . 

1=1.12) 

PUNCH  103 

PUNCH  101. {Al  I  ) ,1 

=1,12) 

DC  73  K  =  13  ,] 7 

73  SI (K)=CI 

DC  75  1  =  1.12 

J«30-l 

75  SI ( J)  =  c  I  (  I  ) 

=l./28. 

BNM 1  )  =  1  .0 

DC  76  1=2,29 

J  =  I-1 

II)=BM(J)-DI> 

DC  77  1=1.29 

7  7  BM<  I  )=bl-i(  1  )/si  (  I  ) 

PUNCH  1     <* 

P'JNCH  101.  (E  'ill, 

1=1 ,29) 

AA=0.0 

DC  78  I=2»28.2 

78  AA*AA+BM(  I  ) 

BB  =  -  • 

DC  79  1=3,27.2 

7  9  BB=BB+BM< I ) 

AA*(BM<29)+4.*AA+2.*BB  )*DX/3. 

X  =  0 

DC  80  1=1,29 

T( I |=BM< I >*X 

P.  X=X+DX 

BB=  . 

DC  81  1=2,28,2 

81   B=BB+T(  I  ) 

cc=o. 

DC  82  1=3,27,2 

82  CC=Cr+T( I ) 

ITC11+T (29)+4. 

• 2.*CC )*DX/3. 

TB=BB/14U. 

TA=AA-TB 

CF=TB/TA 

ST=TA/(TA**2-Tfi**2 ) 

PUNCH  105 

PUNCH  106.  CF.ST 

DC  83  1=1,27,2 

J=  I  +  1 

■ '  -1+2 

T ( I )=<5.*BM( I )+S. 

*BM(J)-BM(K) )*DX/12. 
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A(1)=0.I 

DC  84  1=2,29 

J=  1-1 
84  A ( I )=A( J) +T( J) 

PUNCH  ]■  / 

PUNCH  101, (A( 1 ) ,1  =  1.29) 

DC  8  5  I =1,27,2 

J=I  +  1 

K=l  +7 

T  (  I  )  =  (  5.*A(  I  |+8.*A  (  J  )-A((C)  )*DX/12. 
8  5  T ( J )=( 5.*A(K )+8.*A ( J )-A( I ) |*DX/12. 

B(1)=0.0 

DC  86  1=2,29 

J  =  I-1 
8  6  d(  I  )=i;  (  Jl+Tl  J) 

PUNCH  10  8 

PUNCH  1  1  ,  (l (  i  ) ,1  =  1 ,29) 

DM=B(29)/28. 

DY  =  0 

DC  87  1=1,29 

B ( I )=DY-8( I ) 
8  7  DY=DY+DM 

DC  88  1=1,29 

J=30-I 
8  8  A( J)=B( 1 ) 

DO  89  1=1,29 
89  711 >=ST*(B( I )-CF*A( I  )  ) 

PUNCH  V   9 

PUNCH  1  1  , (T( I  )  ,1  =  1,29) 

100  FCRMATC/15X»14HHCMENT  INERTIA) 

101  FCRKAT(10X,F20.8) 

102  FORMAT! /16X»12HESSENTRICITY) 

103  FCRMATC/.16X.12HSECTICN  AREA) 
i  i,    FORMAT!  /?  X»3HM/1  ) 

10  5  F0R.v,AT(/3X,17riCARRY-CVER-F  ACTOR,  7X»  91 1ST  1FFNESS) 
V   ft  FCRMAT(2F20.8) 

107  FORMAT I/11X.1 /HTHLTA ( CANT  I LEVEL )  ) 

108  FORMAT  </UX»18HDEFLTN(  CANT  1  LEVEL! ) 

109  FORMAT </10X»18HFEM  INFLUENCE  LINE) 
STOP 

END 
DAT* 
in.i  '  .    2778    r  .  10.     7<   .833 
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3UTPUT 

C   C   CALCULATION  OF  HUNCH   ;  I     r;ut  TC  CABLE  STIFFENING     JIN-SIEI 

■-"F«T  INFRTIA 
299*9235' 
222*6649  0 
161*4442  . 
116.2617  i 

67.1  !  71  541  0 

66.740131  I 

50.7460,: 

39.05945700 

31.6804241  0 

26.' 

2  3  .  5  43  2  4« 

21.51056 

FSSFNTRICITY 

1 »39545270 

.96S06436 
175587025 
.54267611 
.45636405 
.35717016 
.25795627 
.1 5874237 
.i~i 96085 
.  B130727 
.u4065369 

ShCI ISN  AREA 

1  .  /  ■ 

.  .  /5i  06ou0 
10*75    01 
10i7500( ' 

i  c .  7 u 

10*50  41 

I     15.1  04001 
!    '  .  r  ' 

10*5'  041     ' 

1     125  ■  •' 
10*25 002 <     0 
.     - 

M/I 

.u^33541fe 
k 00433066 

.    /■};,■■ 

.  '767972 
98  1897 
.'.  1  230787 
.  1' 48327 
.  1920151 

..:/^466r' 

.  :".-i'v66-} 

.02730538 

.02822534 


.1  274 

.  257147] 

.024 

. 

57i77 
.  1826346 

.'■■1. .16966 

.0]  193525 

'.1864 

.00422271 
.00267563 

53983 
92157 

16040 
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CARRY-CVER-FACTCR 

.73432514 


STIFFNESS 

2.03819740 


I'rltTAfCAMILLVEL) 
.     ■.. 

.  '898522 
.06237516 

.    [»';7?f 

.14675515 
,20182786 

.26278166 
.  I  >649l  5 
.42296633 
.54254993 
.62283828 
.76238021 
,841  94331 
.9738'  26? 
1.  ir ! 7120 
i  .16088730 
1.22654]  ' 
1  .32395680 
l  ,37978! 
1.4474  1 ! 
1.49  26  ! 
1.53104400 
1 .56886010 
1 .58591290 
1 .62.  36630 
1.6266  7  i 
l  .65970510 
1.661161* 
1.69377361 


OEFLTMC 


ANTILEVEL ! 

,03729443 

.23052676 

.61885126 

1.21B12570 


-  .  B  7 1 32 7 0 

3.24620610 

4.78292880 

6.72  113 

9.150277  0 

I  "  .  L  L  2 1 1  0 

15.56857! 

19.60229100 

24.16477700 

29.23783200 

34.77383700 

40.7632< 

47.1  56841 

53.93352300 

61  .'  1  ■  «47r  ^ 

68.^7    397    I 

75.92992600 

83.67995600 

91.55967900 

99*56811800 

1l7.67457000 

115.87937000 

124.16855 

I32«.54290OO0 

FEM    INFLUENCE    I  [NE 

.     •         449( 

4.12319140 

8.05597950 
11.71 742400 
15.09609500 
18.0717.4'  [ 
20.64875600 
22.1 

24.12055200 
24.88541000 
25.054242O0 
24.53447 
23.*?67«  f 
21 .94&19300 
2  .  c^62200 
17.91036600 
15.6O71U500 
13.27153500 
10.97120500 

8.89418800 

6.98979010 

5.41564740 

4.0^7611  1  ■ 

2.99044530 

2.10006480 

1 .41 36102C 
.84244487 
.3  -1292  53 

-.00006115 
5TCP   END  CF  PROGRAM  AT   STaTFJ  TOT   ] 09  +  Ol  LINES 
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A  cable  stiffened  structure  is  a  recent  innovation  in 
...  architectural  structure  design.   It  has  the  following 
advantages  over  common  type  of  structural  design:   (1)  re- 
ducing the  dead  load  of  tho  structure,  (2)  efficient  utili- 
zation of  materials  and  (3)  larger  span  or  space  can  be  built. 
The  cable  stiffened  structures  generally  are  high  degree 
inderr.inate  structures.   There  are  many  possible  approaches 
to  the  analysis  of  this  type  of  structure.   In  this  study, 
a  typical  long  span  rigid  frame  with  cable  stiffening  is 
chosen  to  illustrate  some  method  of  stress  analysis.   Two 
possible  methods  which  have  never  been  applied  to  the  cable 
stiffened  structural  analysis  are  introduced  in  this  analysis. 

The  first  method  is  the  application  of  the  principle 
of  superposition.   It  is  based  on  Castigliano' s  theorem  that 
the  deflection  at  each  cable  attachment  can  be  superposed 
by  the  deflections  due  to  external  load  and  due  to  cable  stress. 
So  a  linear  equation  for  deflections  with  functions  of  cable 
stress  can  be  formed  for  each  cable  attachment.   The  number 
of  cables  it  has  is  equal  to  the  number  of  equations  which 
can  be  formed.   Therefore,  cable  stress  can  be  found  by  solv- 
ing the  simultaneous  equations,  the  digital  computer  can  be 
utilized. 

The  second  method  is  applying  the  haunch  member  theory 
to  the  cable  stiffened  structural  analysis.   Because  the  cable 
stiffened  horizontal  girder  of  the  long  span  rigid  frame  can 
be  regarded  as  a  haunched  beam.   The  girder  is  divided  into 


many  small  equal  interval  sections.   In  combining  the  girder 
section  with  the  transformed  section  of  cables,  the  moment  of 
inertia  for  each  section  of  the  girder  can  be  found  and  a 
frame  with  variable  moment  of  inertia  is  formed.   Therefore, 
the  cable  stiffened  rigid  frame  can  be  analyzed  as  a  haunched 
frame . 

Numerical  examples  are  given  to  illustrate  the  methods. 
There  are  examples  with  same  dimension  and  loading  conditions 
for  both  methods.   A  comparison  of  the  results  of  the  two 
methods  shows  the  result  from  the  haunched  frame  is  about  8% 
smaller  than  the  result  worked  by  the  method  of  superposition. 
This  is  because  of  the  haunch  member  is  calculated  section  by 
section,  error  is  introduced  from  these  smaller  sections. 
The  smaller  the  section  is  divided,  the  more  accurate  result 
is  given. 

Larger  space  and  light  weight  of  construction  are  needed 
in  modern  architectural  design.   Cable  stiffening  enables  this 
kind  of  structure  to  be  built.   The  cable  stiffened  rigid 
frame  can  be  designed  for  many  purposes,  such  as  the  construc- 
tion of  large  spaced  field  house,  auditorium,  factory  and 
hangar . 


